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ABSTRACT 


The  theoretical  equilibria  resulting  from  corrosion  of  re¬ 
fractory  nozzle  materials  by  hot  combustion  gases  have  been  determined 
for  a  number  of  systems  of  interest.  The  refractories  considered  have 
include  graphite  (C) ,  tungsten  (W) ,  TaC,  HfN,  TiC,  ZrC,  ZrB2»  MgO,  and 
Hf(>2.  The  bulk  of  the  corrosion  studies  have  been  performed  for  "pure" 
species  of  arbitrary  composition,  including  AIF3,  BF2,  BF3,  BOF,  HB02, 

BeF2»  CO,  CO2,  HC1,  HF ,  H2 ,  H2O,  LiF,  N2  end  condensed  AI2O3  and  LeO. 

Several  real  propellant  systems  also  have  been  investigated.  In  gen¬ 
eral,  corrosion  effects  have  been  determined  for  temperatures  in  the 
range  1500-5000°K,  and  for  a  nominal  total  pressure  (usually  1000  psia) . 

The  results  of  these  studies  have  been  interpreted  in  terms  of 
a  "saturation  number",  defined  as  the  mass  of  wall  material  required  to 
saturate  a  unit  mass  of  reacting  gases.  Saturation  numbers  have  been 
determined  both  for  "Isothermal"  and  for  "adiabatic"  conditions.  Under 
1  certain  Idealized  conditions  the  saturation  number  divided  by  the  density 

of  the  solid  is  directly  proportional  to  the  linear  rate  of  regression  of  the 
|  surface  under  attack,  and  the  "adiabatic  saturation  temperature"  becomes 

1  the  steady-state  temperature  of  the  surface,  allowing  for  corrosion  heat 

effects. 

'  Comparing  the  different  refractories  on  this  basis,  tungsten  (W) 

‘  and  graphite  (C)  generally  rank  the  most  resistant  of  those  materials  con¬ 

sidered.  Tungsten  has  an  advantage  over  graphite  for  combustion  systems 
I  containing  oxygen  or  hydrogen,  while  the  reverse  is  true  for  systems  based 

upon  fluorine.  Tantalum  carbide  (TaC)  also  is  reasonably  resistant  to  at¬ 
tack  by  many  systems.  No  material  considered,  however,  was  capable  of 
withstanding  high-temperature  attack  by  severely  oxidizing  atmospheres, 
such  as  CO2  and  H2O. 

A  limited  comparison  between  measured  rates  of  graphite  nozzle 
erosion  and  rates  predicted  from  a  simplified  equilibrium  model  also  has 
been  made.  At  the  wall  temperatures  typical  of  short  firings  with  state- 
of-the-art  solid  propellants,  reaction  rates  apparently  are  not  sufficiently 
rapid  to  achieve  equilibrium;  the  predicted  erosion  rates  were  high  by  an 
order  of  magnitude  or  so.  The  simple  equilibrium  model  may  show  somewhat 
greater  promise  for  higher-temperature  systems,  but  even  in  this  case  it  no 
doubt  would  have  to  be  modified  to  allow  for  non-unity  Lewis  and  Prandtl 
numbers . 


Theoretical  performance  calculations  previously  conducted  for 
unclassified  and  classified  liquid  bi-propellant  systems  also  are  cited 
and  referenced. 
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SECTION  1 

INTRODUCTION 


This  Is  Volume  I  of  the  Final  Report  In  a  program  which  has 
had  two  objectives,  the  first  being  the  theoretical  appraisal  of  various 
refractory  rocket  liner  materials  from  a  thermochemical  standpoint,  and 
the  second  being  the  computation  of  theoretical  performance  of  certain 
rocket  propellant  combinations.  Chemical  corrosion  studies  reported  In 
the  First*,  Second*,  Third  ,  Fourth  ,  and  Fifth3  Quarterly  Reports  are 
summarized  in  this  volume.  Reference  here  also  is  given  to  the  propellant 
performance  studies  presented  in  Volumes  11^  and  III7  of  the  Second  Quar¬ 
terly  Report.  The  thermodynamic  data  developed  for  new  heavy  metal  com- 

12  3 

bustion  species  and  presented  In  the  first  three  Quarterly  Reports, ’  ’ 
on  the  other  hand,  represent  a  possible  area  of  separate  interest,  and 
accordingly  have  been  compiled  In  Volume  II  of  this  report. 

In  Section  2  of  this  volume,  the  assumptions  and  theory  leading 
to  the  simplified  equilibrium  corrosion  model  are  reviewed,  including 
derivation  of  the  "saturation  number"  as  an  approximate  measure  of 
potential  chemical  attack.  The  computer  program  employed  at  Aeronutronic 
to  calculate  such  numbers  for  both  "Isothermal"  and  "adiabatic"  conditions 
is  described,  and  sources  of  thermodynamic  data  (apart  from  those  compiled 
in  Volume  II)  are  given.  The  individual  corrosion  systems  investigated 
during  the  course  of  this  work  are  reviewed,  Indicating  those  refractory 
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materials  considered  with  certain  "pure"  species  and  those  considered 
with  real  propellant  systems. 

In  a  departure  from  previous  presentations  of  theoretical 
corrosion  results,  different  refractory  materials  are  compared  for  given 
corrosive  atmospheres  in  Section  3  of  this  volume.  In  general,  nine 
separate  refractories  are  considered  in  combination  with  up  to  fourteen 
different  "pure"  species,  from  which  results  relative  ratings  are  derived 
and  discussed.  Comparisons  also  are  shown  for  rocket  nozzle  materials 
attacked  by  five  real  propellant  combustion  systems. 

The  comparisons  between  rates  of  nozzle  erosion  predicted  by 

a 

the  simple  equilibrium  model,  calculated  by  Jones  and  Delaney  ,  and 
measured  by  NASA-Lewis^,  previously  presented  in  the  Fifth  Quarterly 
Report5,  are  reviewed  in  Section  4  of  this  volume.  Theoretical  per¬ 
formance  calculations  also  previously  performed  are  cited  in  Section  5. 
Errata  for  previous  reports  are  given  in  the  Appendix. 
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SECTION  2 

CHEMICAL  CORROSION  STUDIES 


2.1  EQUILIBRIUM  SATURATION  MODEL 

2.1.1  GENERAL  CONSIDERATIONS 

The  theoretical  studies  of  chemical  corrosion  of  refractory 
rocket  liner  materials  by  hot  combustion  gases,  as  carried  out  under 
this  contract,  have  been  based  upon  the  assumption  of  complete  thermo¬ 
chemical  equilibrium  at  the  gas-solid  Interface.  Kinetically  limited 
reactions,  such  as  are  known  to  occur  at  lower  temperatures,  both  homo¬ 
geneously  In  the  gas  phase  and  heterogeneously  with  the  refractory  wall, 
thus  have  been  excluded  from  the  analysis.  Physical  erosion  by  virtue  of 
such  effects  as  aerodynamic  shear,  thermal  shock,  and  liquid  run-off  like¬ 
wise  has  not  been  considered.  While  it  is  recognized  that  the  simplified 
model  assumed  may  lead  to  errors  of  as  much  as  an  order  of  magnitude  or 
more  in  certain  situations,  it  is  felt  that  this  approach  may  be  of  value 
as  a  preliminary,  relative  measure  of  corrosion  resistance  for  use  in 
screening  purposes. 

Under  the  assumption  of  thermodynamic  equilibrium,  the  stability 
of  a  refractory  material  exposed  to  a  given  combustion  atmosphere  may  be 
expressed  in  terms  of  the  quantity  of  refractory  required  to  "saturate"  a 
unit  mass  of  reacting  combustion  products.  The  term  "saturation"  is  meant 
to  encompass  all  possible  gas-solid  or  condensed  phase-solid  reactions 
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forming  either  gaseous  or  condensed  reaction  products,  plus  any  simple 
vaporization/dissociation  of  the  material  in  question.  Indeed,  if  the 
overall  "vapor  pressure"  of  the  refraction  material  by  itself  should 
exceed,  at  a  given  temperature,  the  total  pressure  of  the  system,  satura¬ 
tion  could  never  be  reached.  (Thus,  unstable  ablative  plastics  are  not 
covered  by  the  model.)  On  the  other  hand,  it  is  possible  for  a  combustion 
mixture  containing  the  same  elements  as  the  refractory  to  exceed  satura¬ 
tion  in  advance,  so  that  deposition  rather  than  corrosion  theoretically 
would  be  expected.  Just  as  negative  erosion  would  be  predicted  for  the 
case  just  cited,  positive  or  net  erosion  would  be  calculated  for  corrosion 
reactions  forming  solid  products,  even  though  such  products  in  reality 
might  form  a  protective  layer  over  the  refractory  surface. 

2.1.2  "ADIABATIC"  SATURATION  CONDITIONS 

The "saturation  number"  of  a  refractory  material  in  a  given  mix¬ 
ture  (i.e.,  the  mass  of  wall  material  required  to  saturate  a  unit  mass  of 
combustion  products)  can  be  defined  either  for  "isothermal"  or  for  "adiabatic" 
conditions,  and  has  been  denoted  by  "A"  and  by  "B" ,  respectively,  in  previous 
reporting  under  this  contract.  More  generally,  however,  the  designation  "B" 
is  used  for  the  saturation  number  under  conditions  not  necessarily  specified. 
The  definition  of  the  "isothermal"  saturation  number  "A"  for  a  fixed  tem¬ 
perature  and  pressure  is  unequivocal,  depending  only  on  basic  thermo¬ 
dynamic  considerations.  The  concept  of  the  "adiabatic"  saturation  number, 
on  the  other  hand,  does  involve  a  relationship  between  bulk  flow  conditions, 
saturation  conditions  at  the  wall,  and  rates  of  heat  and  mass  transfer  within 
the  boundary  layer  for  a  given  geometry,  and  hence  depends  on  additional 
assumptions  of  somewhat  more  questionable  validity.  A  simplified  derivation 

of  equations  treating  the  "adiabatic  saturation"  case  was  presented  in  the 

2 

Second  Quarterly  Report  ,  and  will  be  reviewed  here. 

Under  the  assumptions  of  unity  Lewis  and  Prandtl  numbers, 
as  considered  by  Denison^®’  ^  and  by  Bartlett^.,  and  under 
certain  other  conditions,  including  equilibrium  at  the  wall  with 
blowing  and  net  heat  flow  neglected,  the  wall  temperature  in  an 
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adiabatic  equilibrium  flow  system  at  constant  pressure  becomes  equal  to  a 
temperature  commonly  designated  as  the  "adiabatic  saturation"  temperature, 
and  defined  by  an  enthalpy  balance  between  free  stream  (stagnation)  gas, 
virgin  wall  material,  and  saturated  mixture  adjacent  to  the  wall.  Depending 
on  the  heat  effect  of  the  reaction  involved,  be  it  oxidation  or  simple 
"vaporization",  the  adiabatic  saturation  temperature  may  be  either  higher 
or  lower  than  that  corresponding  to  free  stream  stagnation.  The  balance 
to  determine  the  adiabatic  wall  conditions  may  be  written  as  follows: 

(1  +  B)h  -  h  +  Bh  (1) 

'  sw  se  w 

in  which 

B  -  adiabatic  saturation  number,  gms  wall  material  to 
saturate  one  gm  of  free  stream  gas  at  Tsw,  the 
saturated  wall  temperature. 

hflw  -  enthalpy  of  saturated  mixture  at  T8W,  kcal/100  gms. 
hse  »  enthalpy  of  free  stream  gas  at  Tse,  the  stagnation 
temperature,  kcal/100  gms. 

1^  -  enthalpy  of  virgin  wall  material  at  the  temperature, 

Tw,  at  which  it  is  added  to  the  Bystem,  kcal/100  gms. 

2.1.3  ACTUAL  RATES  OF  EROSION 

Whereas  the  determination  of  the  adiabatic  saturation  condition 

from  equation  (1)  depends  only  on  relative  rates  of  heat  transfer  versus 

mass  transfer,  absolute  transport  rates  must  be  estimated  in  order  to 

predict  an  actual  rate  of  surface  erosion.  Again  considering  the 

assumptions  of  unit  Lewis  and  Prandtl  numbers  in  the  manner  of  Denison^0’!!, 

*  2 

it  was  shown  in  the  Second  Quarterly  Report  that  the  overall  erosion  rate 
at  steady  state  should  be  proportional,  for  gaseous  reactants,  to  the 
adiabatic  saturation  number  "B": 

*  “  <3F)b  <2> 

*The  model  described  here,  and  designated  as  the  "simplified  Denison"  approach 
should  not  be  confused  with  the  much  more  general  case  originally  developed 
by  Denison. 
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in  which 

2 

m  -  erosion  rate,  lb/hr-ft  . 
hc  ■  overall  heat  transfer  coefficient,  Btu/hr-ft2-°F. 
Cp  ■  gas  heat  capacity,  Btu/lb-°F. 

The  linear  regression  rate  "r"  then  follows  directly: 


in  which 

t  -  regression  rate,  ft /hr. 

»  density  of  wall  material,  lb/ft  . 

In  a  situation  with  net  heat  flow  into  the  refractory  wall,  the  balance 

represented  by  equation  (1)  would  no  longer  hold  true.  Nevertheless,  it 

can  be  shown  that  equations  (2)  and  (3)  should  still  apply,  provided  that 

"B"  is  understood  in  this  case  to  be  merely  the  saturation  number  at  some 

wall  temperature  T#w*  which  must  be  known  before  the  erosion  rate  can  be 

estimated.  This  approach  was  adopted  in  comparing  predicted  versus  measured 

erosion  rates  in  the  Fifth  Quarterly  Report5  under  this  contract. 

Assuming  the  approximate  validity  of  equations  (2)  and  (3)  for 

Lewis  numbers  near  unity,  perhaps  the  only  remaining  controversial  aspect 

of  calculations  to  predict  erosion  rates  lies  in  the  definitions  of  "h" 

he 

and  "Cp",  or  more  precisely,  of  their  ratio  (r-)  .  For  non-dissociating, 

Lp 

non-reacting  systems,  the  overall  heat  transfer  coefficient  "hc"  is  un¬ 
ambiguous  in  meaning,  and  may  be  estimated  using,  for  example,  the  Bartz^ 
equation  for  rocket  nozzle  heat  transfer.  For  the  dissociating  case,  the 
approach  taken  here  in  rough  calculations  has  been  to  assume  that  both 

"hc"  and  "Cp"  increase  in  about  the  same  ratio,  so  that  the  enthalpy  heat 

he 

transfer  coefficient  (— )  may  still  be  estimated  from  standard  correlations. 

up 

While  more  sophisticated  heat  transfer  analyses  are  available,  their  use 
with  the  simple  equilibrium  saturation  model  has  not  been  considered 
justified  in  view  of  other  approximations  required. 
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Aside  from  Che  question  of  attainment  of  equilibrium  itself,  the 
assumption  of  unity  Lewis  numbers  in  deriving  equations  (1),  (2)  and  (3) 
has  been  considered  the  most  likely  source  of  error,  A  preliminary  and 
elementary  treatment,  applying  corrections  to  these  equations  for  blowing, 
heat  flow  into  the  wall,  and  non-unity  Lewis  numbers  has  been  under 
development^,  but  has  not  been  tested  against  experimental  data.  In  view 

of  the  uncertainties  present  in  this  model,  it  is  not  presented  here. 

2 

In  the  Second  Quarterly  Report  it  was  shown  that  the  saturation 
number  "B"  for  a  combustion  mixture  without  condensed  species  would  be 
equivalent  to  the  "blowing  parameter"  (B1)  defined  by  DenlsonlO,  if,  in 
addition  to  unity  Lewis  and  Prandtl  numbers,  the  Reynolds  analogy  between 
heat  and  momentum  transfer  was  assumed.  In  other  words: 


in  which 


B  *  B' 


e.vf> 


(A) 


B'  -  blowing  parameter,  dimensionless. 

■  gas  density  at  outer  edge  of  boundary  layer,  lb/ft  . 
Ug  -  gas  velocity  at  outer  edge  of  boundary  layer,  ft /hr. 
f  -  Fanning  friction  factor,  dimensionless. 


2.2  COMPUTER  PROGRAM 

2.2.1  NOZZLE  ABLATION  PROGRAM 

The  computation  of  isothermal  multicomponent  equilibrium  with 
an  excess  of  a  given  condensed  phase  may  be  accomplished  with  a  number  of 
available  programs,  including  the  single -point  option  of  the  Aeronutronlc 
Rocket  Propellant  Performance  Program  written  for  the  IBM  7090  computer. 

More  conveniently,  however,  a  separate  equilibrium  program  might  be  used 
which  would  automatically  converge  to  a  saturated  composition,  and  then 
compute  the  value  of  "A",  the  isothermal  saturation  number.  Such  a  program, 
known  as  the  Nozzle  Ablation  Program,  was  written  for  the  Aeronutronlc 
IBM  7090  computer  and  was  described  in  the  First  Quarterly  Report*  under 
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this  contract.  The  procedure  used  was  simply  to  add  2000  gms  of  wall 
material  to  100  gins  of  the  system  under  consideration,  and  then  totompute 
the  resulting  equilibrium  state,  assuming  that  saturation  would  be  achieved. 

With  the  advent  of  the  adiabatic  saturation  model  for  describing 

steady-state  chemical  corrosion  of  rocket  liner  materials,  the  original 

Nozzle  Ablation  Program  had  to  be  rewritten  to  Include  iterative  solution 

of  the  enthalpy  balance  represented  by  equation  (1).  As  described  in  the 

2 

Second  Quarterly  Report  ,  convergence  is  approached  by  assuming  trial 
temperatures,  computing  isothermal  saturation  conditions  for  each  tem¬ 
perature,  and  then  changing  temperature  based  on  the  difference  between 
right  and  left  hand  sides  of  equation  (1).  The  program  was  also  modified 
to  add  200,  then  400,  then  800,  and  finally  1600  gms  of  wall  material  to 
100  gms  system  until  arriving  at  isothermal  saturation,  in  order  to  avoid 
the  large  excess  of  refractory  phase  often  present  in  the  final  Bystem  if 
a  fixed  addition  ratio  of  20:1  was  used. 

The  program  as  presently  written  has  several  options  with  regard 
to  choice  of  free  stream  or  saturated  wall  temperatures.  For  a  real  propel¬ 
lant  system  with  overall  heat  of  formation  and  throat  pressure  known  in 
advance,  the  Nozzle  Ablation  Program  can  compute  the  throat  stagnation  con¬ 
dition  in  exactly  the  same  way  that  a  chamber  condition  is  computed  by  the 
Rocket  Propellant  Performance  Program,  and  then  can  calculate  both  iso¬ 
thermal  and  adiabatic  saturation  numbers  for  that  "free  stream"  tempera¬ 
ture.  In  the  option  normally  used  in  calculations  for  atmospheres  of 
arbitrary  composition  under  this  contract,  however,  free  stream  temper¬ 
atures  in  500°K  Increments  from  1500°K  up  to  5000°K  are  assumed,  and  then 
both  an  isothermal  "A"  and  an  adiabatic  "B"  are  computed  for  eech  such 
temperature .  (The  computations  thus  are  arbitrarily  extended  without 
modification  to  temperatures  above  the  melting  point  of  the  refractory  in 
question.)  Two  options  also  are  available  for  choosing  the  initial  wall 
temperature  "Tw"  used  to  define  "h^',  the  wall  material  enthalpy  entering 
into  equation  (1):  (a)Tw  “  298°K,  (b)Tw  »  T8W,  the  adiabatic  saturated 
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wall  temperature.  As  a  possibly  more  realistic,  and  certainly  more  severe 
assumption,  the  second  case  has  been  adopted  in  all  calculations  reported 
under  this  contract. 

The  usual  assumptions  common  to  most  thermochemical  equilibrium 
programs,  i.e.,  ideal  behavior  in  the  gas  phase  and  no  solid  or  liquid 
solutions  formed  between  condensed  phases,  also  have  been  incorporated  into 
the  Nozzle  Ablation  Program.  As  a  result,  the  program  will  not  converge 
satisfactorily  for  systems,  particularly  occurring  with  condensed  reactants, 
which  form  only  condensed  products  at  the  pressure  under  consideration. 

In  some  cases  this  difficulty  can  be  avoided  by  the  inclusion  in  the  input 
of  small  amounts  of  some  inert  gas  such  as  argon,  but  even  then,  if  species 
derived  from  the  wall  material  are  not  present  in  the  gas  phase  at  equil¬ 
ibrium  in  significant  amounts,  convergence  will  not  be  achieved.  Such 
cases  normally  indicate  stoichiometric  reactions  which  can  be  best  treated 
by  simple  hand  calculations. 

2.2.2  ESTIMATION  METHOD 

In  order  to  permit  convenient  estimation  of  saturation  numbers 
for  refractories  exposed  to  complex  combustion  systems,  a  method  was  de¬ 
sired  by  which  these  numbers  might  be  derived  empirically  from  the  satura¬ 
tion  numbers  determined  under  this  contract  for  individual  "pure"  species, 
without  need  to  perform  separate  computer  calculations  for  each  new  case. 
Such  a  method,  based  upon  the  concept  of  "additivities",  was  Investigated 

4 

in  the  Fourth  Quarterly  Report  .  It  was  shown  that  the  Isothermal  satura¬ 
tion  numbers  "A"  were  additive  on  a  weight  basis  for  those  corrosion  re¬ 
actions  which  (a)  involved  no  net  change  in  gaseous  moles,  (b)  consisted 
of  simple  vaporlzation/decomposition  of  the  refractory,  or  (c)  went  to 
stoichiometric  completion.  For  individual  reactions  involving  either  a 
net  increase  or  a  net  decrease  in  gaseous  moles,  moreover,  algebraic 
equations  could  be  derived  which  corrected  for  the  effects  of  pressure  and 
dilution,  yielding  "blending"  saturation  numbers  which  also  were  additive. 
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The  analysis  did  not  extend  to  the  prediction  of  adiabatic  aaturation 
numbers,  but  vaa  concerned  only  with  the  Isothermal  "A's” . 

In  testing  the  estimation  method  for  true  propellant  systema,  it 
was  found,  however,  that  interactions  occurred  in  many  situations  of  in* 
terest,  such  that  any  prediction  based  upon  additivities  would  necessarily 
fail.  Such  an  interaction  could  lead  either  to  a  decrease  in  saturation 
number  for  the  mixture,  or,  alternatively,  to  an  increase.  As  an  example 
of  the  former,  attempts  to  combine  saturation  numbers  for  graphite  +  H2 
and  graphite  +  F2  would  overpredict  the  corrosivity  of  a  nominal  mixture 
of  H2  +  F2  at  the  same  temperature,  owing  to  failure  to  take  into  account 
the  stability  of  the  gaseous  HF  molecule.  As  an  example  of  the  latter4, 
combining  the  saturation  numbers  for  graphite  +  H2  and  graphite  +  N2  tends 
to  underpredict  the  corrosivity  of  a  mixture  containing  H2  +  N2,  since 
gaseous  HCN  becomes  a  major  reaction  product  of  graphite  with  the  mixture. 
Therefore,  the  estimation  method  has  been  deemed  generally  unsatisfactory 
for  predicting  theoretical  corrosion  with  real  propellant  systems. 

2.3  THERMODYNAMIC  DATA 

The  thermodynamic  data  for  combustion  species  used  in  equilibrium 
computations  under  this  contract  were  taken,  wherever  available,  from  the 
JANAF  tables ^  In  certain  cases,  particularly  for  condensed  species, 
the  JANAF  data  did  not  cover  the  entire  range  from  298°K  up  to  6000°K. 

In  such  instances,  the  JANAF  data  usually  were  extended,  using  best 
estimates  of  heat  capacities,  and,  where  applicable,  heats  of  fusion  and 
melting  points.  In  a  few  cases,  where  the  species  involved  was  known  to 
become  increasingly  unstable  at  higher  temperatures,  the  data  were  trun¬ 
cated  rather  than  extrapolated,  so  that  the  species  would  not  be  considered 
above  a  certain  temperature. 

For  several  heavy  metal  compounds  considered  as  refractories, 
the  pertinent  species  families  either  were  not  adequately  covered  in  the 
JANAF  tables,  or,  in  a  few  cases,  were  not  even  given  at  all  (Ex.  -  Ta,  Hf, 
Mo,  Nb,  Th) .  Thus,  a  significant  portion  of  the  effort  under  this  contract 
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was  devoted  to  compiling  the  thermodynamic  data  necessary  to  consider  some 

of  these  chemical  systems.  Where  available,  data  published  in  sources 

aside  from  JANAF  were  selected,  but  in  most  cases  estimation  was  required, 

using  methods  standard  throughout  the  literature.  Compilations  of  such 

12  3 

data  were  given  in  the  First  ,  Second  and  Third  Quarterly  Reports.  In 

order  to  make  these  data  more  easily  retrievable,  however,  lue  thermodynamic 

properties  of  57  species  of  Ta,  W  and  Hf  not  contained  in  the  JANAF  tables 

have  been  reprinted  as  Volume  II  of  this  report. 

In  addition  to  those  species  listed  in  Volume  II,  the  thermo- 

dynamic  data  previously  estimated  for  WCl, (g) ,WC1, (g) ,WF(g) ,  and  ZrB9fc)  have 

5 15  b  1 

recently  been  superseded  by  new  JANAF  data  .  The  thermodynamic  properties 
of  condensed  HfCl^Cc),  HfCl^(c),  HfF^c),  and  HfF^(c)  moreover  were  taken, 
in  calculations,  to  be  identical,  except  for  molecular  weight,  to  the 
JANAF  data  for  the  corresponding  Zr  species.  The  data  for  five  additional 
species,  C  H(g),  ZrOCl(g),  ZrOCl  (g),  ZrOF(g),  and  ZrOF  (g) ,  were  previously 

*  ^  1£  I?  * 

published  by  Aeronutronic  under  another  contract.  ’  In  order  that  the  exact 
thermodynamic  values  assumed  for  the  chemical  corrosion  studies  be  properly 
available,  the  data  for  C2H(g),  the  three  W  species,  and  the  five  Zr  species 
just  cited  are  reprinted  in  the  Appendix  of  this  volume. 

A  listing  of  all  pertinent  species  considered  in  the  chemical 
corrosion  studies  and  the  propellant  performance  studies  carried  out  under 
this  contract  is  given  in  Table  I.  The  designation  "/c"  after  a  species 
indicates  a  condensed  phase,  and  is  meant  to  include  both  liquid  and  solid 
states,  where  applicable.  (In  the  JANAF  tables,  the  liquid  and  solid  are 
considered  as  separate  species.) 
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TABLE  I 


IN  EQUILIBRIUM  CALCULATIONS 


C  Specie* 

C  Specie* 

(eont.) 

(cont.) 

c5 

ch3f 

CO 

ch2f2 

co2 

chf3 

CH 

HCOF 

ch2 

CNF 

ch3 

CC1 

ch4 

CC14 

c2h 

coci2 

C2H2 

ch3ci 

c2h4 

CNC1 

CHO 

ccif3 

ch2o 

cci2f2 

CN 

COC1F 

C2N2 

chcif2 

c4n2 

C/c 

HCN 

HCNO 

B  Species 

CF 

B 

cf2 

BO 

cf3 

bo2 

cf4 

b2o2 

c2f2 

b2o3 

c2f4 

BH 

cof2 

bh2 

o 
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SPECIES  CONSIDERED  IN  EQUILIBRIUM  CALCULATIONS 


Species 
(cont .) 

A1  Species 

LI  Species 

bh3 

A1 

Li 

HBO 

AlO 

Ll2 

hbo2 

ai2o 

LiO 

h2bo2 

ai2o2 

li2o 

H3BO3 

A1H 

L1H 

BF 

A10H 

LlOH 

bf2 

haio2 

(lioh)2 

BF3 

AIN 

LIN 

bof 

A1F 

L1F 

(bof)3 

aif2 

(lif)2 

hbof2 

AIF3 

(lif)3 

BC1 

A10F 

LlOF 

BC1, 

A1C1 

L1C1 

BCl, 

aici2 

(licd2 

BOC1 

AICI3 

HC 10 

(BOCI)3 

A10C1 

L12C1F 

BCIF 

A1C1F 

libo2 

bcif2 

aicif2 

liaif4 

bci2f 

aici2f 

Li/c 

®2cia 

Al/c 

L120/c 

BC 

A1203/c 

LiH/c 

B/c 

AlN/c 

LlOH/c 

B203/c 

A1F3/c 

LI3N/C 

HB02/c 

AICI3/C 

LIF/c 

H3BO3/C 

A1aC3/c 

LiCl/c 

BN/c 

L12C2/c 

Lt  Species 
(cont.) 

L1B02/c 

L12Ba07/c 

l12b6°10/c 

Li2Bg0^3/c 

LiA102/c 

Be  Species 

Be 

BeO 

(BeO)2 

(BeO) 3 

(BeO)A 

(BeO)j 

(BeO)6 

BeH 

BeH2 

BeOH 

Be (OH) 2 

BeF 

BeF2 

BeCl 

BeCl2 

BeCIF 

BeC2 

BeB20A 
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SPECIES  CONSIDERED  IN  EQUILIBRIUM  CALCULATIONS 


Be  Species 

Mg  Species 

(cont . ) 

(cont . ) 

LiBeF3 

MgF2/c 

Be/c 

MgCl2/c 

BeO/c 

MgC2/c 

Be (OH) 2 /c 

Mg2C3/c 

Be3N2/c 

Al2Mg04/c 

BeF2/c 

BeCl2/c 

Ta  Species 

Be2C/c 

Ta 

Al2Be04/c 

TaO 

TaO- 

Mr  Species 

A 

TaF 

Mg 

TaF2 

MgO 

TaF  3 

MgH 

TaF4 

MgOH 

TaF5 

MgF 

TaOF 

MgF2 

TaOF2 

MgCl 

TaOF3 

MgCl2 

Ta02F 

MgCIF 

TaCl 

Mg/c 

TaCi2 

MgO/c 

TaCl3 

MgH2/c 

TaCl4 

Mg(OH)2/c 

T«C15 

Mg3N2/c 

TaOCl 

TaOCl2 

Ta  Species 

Ti  Species 

(cont .) 

(cont . ) 

TaOCl 3 

Ti203/c 

Ta02Cl 

Ti305/c 

Ta/c 

TiN/c 

Ta205/c 

TiF2/c 

TaN/c 

TiF3/c 

TaCl3/c 

TiF4/c 

TaCl4/c 

TiCl2/c 

TaClj/c 

TiCl3/c 

TaC/c 

TiC /e 

TaB2/c 

TiB2/c 

Li2Ti03/c 

Ti  Species 

MgTiO-j/c 

Ti 

MgTi205/c 

TiO 

Mg2Ti04/c 

Ti02 

TiF 

Zr  Species 

TiF2 

Zr 

TiF3 

ZrO 

TiF4 

Zr02 

TiCl 

ZrF 

TiCl2 

ZrF2 

TiCl3 

ZrF3 

TiCl4 

ZrF4 

Ti/c 

ZrOF 

TiO/c 

ZrOF2 

Ti02/c 

ZrCl 
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SPECIES  CONSIDERED  IN  EQUILIBRIUM  CALCULATIONS 


Zr  Species 

W  Species 

(cont. ) 

(cont.) 

ZrCl2 

h2wo4 

ZrCl3 

WF 

ZrCI4 

wf2 

ZrOCl 

wf4 

ZrOCl2 

wf5 

Zr/c 

wf6 

Zr02/c 

wof2 

ZrN/c 

wof4 

ZrF2/c 

wo2f2 

ZrF3/c 

WC1 

ZrF4/c 

wci2 

ZrCl2/c 

WC14 

ZrCl3/c 

WC15 

ZrCl4/c 

WC16 

ZrC/c 

WOCI2 

ZrB2/c 

woci4 

W02C12 

W  Species 

W/c 

W 

W02/c 

WO 

W03/c 

W02 

H2W04/c 

WO- 

WF6/C 

J 

(W03)3 

WOF4/c 

(wo3)4 

WC12/c 

W  Species 

Hf  Species 

(cont. ) 

(cont.) 

WCl4/c 

Hf/c 

WClj/c 

Hf02/c 

WCl6/c 

HfN/c 

WOCl4/c 

HfF2/c 

W02Cl2/c 

HfF3/c 

WC/c 

HfF4/c 

WB/c 

HfCl2/c 

HfCl3/c 

Hf  Species 

HfCl4/c 

Hf 

HfC/c 

HfO 

HfB2/c 

Hf02 

HfF 

Rare  Gas  Species 

HfF2 

He 

HfF3 

Ar 

HfF, 

Kr 

HfOF 

HfOF2 

HfCl 

HfCl2 

HfCl3 

HfCl4 

HfOCl 

HfOCl2 
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2.4  SYSTEMS  INVESTIGATED 

The  refractory  material*  atudled  during  thla  program  were 
selected  from  an  original  list  Including  the  following*®: 


Metals; 

Carbides; 

Nitrides; 

Oxides: 

Borides; 


C,  W,  Ta,  Mo 

TaC,  HfC ,  NbC,  ZrC,  TIC,  AI4C3,  B4C,  SIC 
TaN,  HfN,  ZrN,  TIN,  AIN,  BN,  S13N4 
MgO,  Zr02,  BeO,  Th02 
HfB2,  ZrB2,  T1B2 


Those  materials  actually  Investigated,  either  through  rigorous  machine 
computation  or  by  approximate  hand  calculation,  were  as  follows:  C,  H, 

TaC,  HfC,  ZrC,  TIC,  HfN,  TIN,  MgO,  Zr02,  Th02,  and  ZrB2. 

The  major  portion  of  the  theoretical  corrosion  studies  were 
carried  out  for  "pure"  combustion  species,  rather  than  for  real  propellant 
systems.  The  species  considered  Included  AIF3,  BF2 ,  BF-j,  BOF,  HB02 ,  BeF2 , 
CO,  C02,  HCl,  HF,  H2,  H20,  LiF,  N2,  A120-j/c,  and  BeO/c.  It  was  recognised, 
however,  that  at  the  "free-atream"  temperatures  and  pressures  Involved, 
many  of  these  nominal  species  would  be  largely  dissociated  or  dlspro- 
portionated. 

The  "pure"  species  systems  actually  investigated  are  shown  in 
Table  II  in  terms  of  a  matrix  Involving  9  refractories  and  16  species. 
Except  as  otherwise  to  be  noted,  all  calculations  were  for  a  total 
pressure  of  1000  psia,  and  for  temperatures  In  the  range  1500-5000°K. 

The  results  were  presented  in  terms  of  two  types  of  plots,  one  type 
showing  isothermal  "A"  and  adiabatic  "B"  as  functions  of  wall  temperature 
and  free-stream  temperature,  respectively,  and  the  other  type  showing  In¬ 
dividual  species  concentrations  in  the  saturated  mixture  as  functions  of 
wall  temperature.  The  relative  corrosivities  versus  temperature  and  the 
heat  effects  involved  could  be  read  from  the  first  type  of  plot,  while 
the  ablation  chemistry  could  be  interpreted  from  the  second  type. 
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TABU  II 

CHEMICAL  CORROSION  OF  REFRACTORY  MATERIALS ; 
"PURE"  SPECIES  SYSTEMS  STUDIED 


C/c 

W/c 

TaC/c 

HfN/c 

TiC/c 

ZrC/c 

ZrB2/c 

MgO/c 

HfO. 

A1F3 

2 

3 

4 

5 

- 

- 

- 

- 

- 

bf2 

2 

3 

4 

5 

- 

- 

- 

- 

- 

bf3 

2 

3 

4 

5 

- 

- 

- 

- 

- 

BOF 

2 

3,4 

4 

5 

- 

- 

- 

- 

- 

hbo2 

2 

(3) 

(4) 

- 

- 

- 

- 

- 

- 

BeF2 

2,6 

3,6 

4,6 

5 

- 

- 

- 

- 

- 

CO 

1,2 

3,4 

4 

5 

2 

2 

2 

- 

- 

co2 

1,2 

3,4 

4 

5 

2 

2 

2 

4 

4 

HC1 

1,2 

3 

4 

5 

2 

2 

2 

- 

- 

HF 

1,2 

3 

4 

5 

2 

2 

2 

- 

- 

h2 

1,2 

6 

4 

5 

2 

2 

2 

- 

- 

h2o 

1,2 

6 

4 

5 

2 

2 

2 

- 

- 

LiF 

2 

3 

4 

5 

- 

- 

- 

- 

- 

N2 

1,2 

6 

4 

5 

2,4 

2,4 

2 

- 

- 

AI2O3/ c 

2,5 

5 

5 

5 

- 

- 

- 

- 

- 

BeO/c 

2 

- 

. 

. 

- 

. 

. 

. 

_ 

KEY:  1  -  First  QPR,  9/15/621  4  -  Fourth  QPR,  6/15/634 

2  -  Second  QPR,  12/15/622  5  -  Fifth  QPR,  9/15/635 

3  -  Third  QPR,  3/15/633  6  -  Final  Report,  12/15/63 
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As  may  be  aeen  from  Table  II,  graphite  was  considered  with 
•even  non-metallized  species  In  the  Firgt}  and  Second^  Quarterly  Reports, 
and  with  the  remaining  metallic  species  also  In  the  Second  Quarterly.  The 
calculations  presented  in  the  First  Quarterly  Report  covered  only  the 
isothermal  saturation  case,  while  In  the  Second  Quarterly,  both  Isothermal 
and  adiabatic  saturation  parameters  were  included.  Also  shown  in  the 
Second  Quarterly  were  the  effects  of  total  pressure  on  theoretical 
corrosion  parameters  for  graphite  exposed  to  H2,  ^0  and  CO2.  Only  in 
the  case  of  the  C/c  +  H2  reaction  was  there  an  appreciable  net  change 
in  gaseous  moles,  and  hence,  a  significant  pressure  effect. 

Theoretical  corrosion  studies  for  tungsten  attacked  by  11  cotmnon 

3 

combustion  species  were  presented  in  the  Third  Quarterly  Report  .  The 
calculations  for  W/c  and  BOF,  CO,  and  CO2  were  revised  in  the  Fourth 
Quarterly  Report^,  based  upon  new  JANAF^  data  for  the  gaseous  tungsten 
oxides.  The  W/c  +  HjO  case  was  again  deferred,  however,  pending  adoption 
of  JANAF  tables  for  gaseous  and  condensed  H2W0^.  It  Is  also  interesting 
to  note  that  whereas  the  system  C/c  +  HB02  yielded  reasonable  saturation 
compositions  over  the  entire  temperature  range^,  the  system  W/c  +  HB02 
apparently  was  completely  condensed  except  at  temper atures  In  the  vicinity 
of  5000°K.^  The  difference  resulted  from  the  calculated  high  stability 
of  the  species  HBO2/C,  which  was  added  to  the  thermodynamic  file  In  the 
interim  between  the  two  reports.  In  view  of  the  considerable  extrapolation 
involved  in  predicting  high  temperature  thermodynamic  values  for  condensed 
HB02/c,  all  calculations  including  this  species  are  considered  extremely 
unreliable . 

The  theoretical  corrosion  parameters  for  TaC/c  attacked  by 
14  "pure"  combustion  species  were  presented  in  the  Fourth  Quarterly 
Report^,  and  corresponding  results  for  HfN/c  considered  as  a  possible  re¬ 
fractory  were  given  in  the  Fifth  Quarterly  Report^.  Again,  the  system 
TaC/c  4  HBO2  was  found  to  give  a  rather  meaningless  result,  and  was  not  re¬ 
ported  in  detail.  Also  for  this  reason,  HB02  as  a  "pure"  species  was  not 
even  calculated  with  HfN/c  as  a  ablator. 
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The  three  refractories  TiC/c,  ZrC/c,  and  Zrl^/c  were  considered 
only  with  the  seven  non-metalllzed  species  CO,  CO2,  HC1,  HP,  H2,  HjO,  and 
N2>  These  data  all  were  presented  In  the  Second  Quarterly  Report?.  In 
the  Fourth  Quarterly  Report^,  the  refractory  oxides  were  considered  briefly 
for  application  in  a  CO2  atmosphere.  Rigorous  theoretical  corrosion  cal¬ 
culations  were  performed  for  MgO/c  and  Hf02/c,  while  rough  hand  calculations 
were  made  for  ThOj/c,  the  highest -melting  oxide. 

Besides  the  normal  chemical  corrosion  calculations  for  1000  psla 
total  pressure,  certain  special  cases  also  were  considered.  In  the  Fourth 
Quarterly  Report^,  hand  calculations  were  employed  in  comparing  the  relative 
stabilities  of  TaC/c,  TiC/c,  and  ZrC/c  in  N2  atmospheres,  in  other  words, 
in  estimating  the  partial  pressures  of  N2  in  equilibrium  simultaneously 
with  the  metal  carbide,  metal  nitride,  and  graphite.  In  the  Fifth  Quarterly 
Report^,  the  results  of  machine  computations  were  give.,  for  the  theoretical 
corrosion  of  C/c,  W/c,  TaC/c,  and  HfN/c  by  condensed  A^O^/c,  under  con¬ 
ditions  such  that  both  condensed  phases  would  be  present  at  equilibrium. 

The  true  total  pressure,  then,  would  vary  with  temperature  in  much  the 
manner  of  a  "mixed  vapor  pressure".  The  single-point  feature  of  the 
Rocket  Propellant  Performance  Program  was  used  for  these  calculations  , 
with  argon  added  to  yield  a  nominal  total  pressure. 

Inasmuch  as  revised  JANAF  data  for  gaseous  and  condensed  H2WO4 
have  recently  been  received?®,  the  W/c  +  H2O  system  was  computed  for  in¬ 
clusion  in  this  report.  The  two  systems  W/c  +  H2  and  W/c  +  N2,  originally 
excluded  owing  to  their  expected  non-reactivity,  also  were  treated,  but 
only  by  simple  hand  calculation,  in  order  to  show  the  effect  of  tungsten 
volatility.  The  species  concentration  plots  for  these  three  cases  are 
presented  in  Figure  1,  in  which  x  (moles  of  gas  per  100  grams  system)  and 

condensed  species  concentrations  (moles  per  100  grams  system)  are  shown 

3 

according  to  the  practice  first  adopted  in  the  Third  Quarterly  Report  . 

The  "A"  and  "B"  saturation  parameters  for  W/c  and  H2,  H2O,  and  N2  are 

g 

presented  (in  terms  of  -75-  )  as  part  of  Figures  12,  13,  and  15,  respectively. 
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FIGURE  1.  SPECIES  PRESENT  UNDER  EQUILIBRIUM  CONDITIONS  AT  1000  PSIA 
FOR  W/c  EXPOSED  TO  H„,  H„0.  AND  N  . 
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As  a  result  of  a  recent  drastic  revision  In  the  JANAF  heat  of 
formation  for  gaseous  BeF,  the  theoretical  corrosion  calculations  pre¬ 
viously  performed  for  C/c,  W/c,  and  TaC/c  (but  not  HfN/c)  exposed  to  BeF 2 
were  Invalidated,  and  also  had  to  be  repeated.  In  Figure  2  are  shown  the 
new  composition  plots  for  these  three  refractories  attacked  by  BeF2‘,  the 
new  "A"  and  "B"  values  (divided  by  densities)  are  indicated  in  Figure  7. 

The  effect  of  the  change  in  thermodynamic  properties  for  BeF  is  to  make 
a  BeF 2  atmosphere  much  less  corrosive  than  was  formerly  believed. 

A  limited  number  of  real  propellant  systems  were  investigated 
for  chemical  corrosion  effects  in  the  course  of  this  program.  In  Section  3 
of  the  Fifth  Quarterly  Report^,  saturation  numbers  were  given  for  graphite 
attacked  by  Arcite  368,  Arcite  373,  "NP",  and  "NF"  propellants.  Calcula¬ 
tions  also  have  been  made  for  W/c  and  TaC/c  with  the  "NF"  propellant,  and 
for  various  refractories  exposed  to  F2/H2  (fuel-rich) ,  "hot  solid",  and 
"BOF  simulation"  mixtures.  The  theoretical  corrosion  parameters  for 
these  latter  systems  are  given  in  Figures  17  through  21  of  this  report. 

In  the  section  to  follow,  refractory  materials  are  rated  by 
their  comparative  resistance  to  various  corrosive  atmospheres.  This  pre¬ 
sentation  is  exactly  the  opposite  of  that  previously  followed  in  the 
Quarterlies,  in  which  different  "pure"  species  were  compared,  at  least 
indirectly,  for  each  separate  refractory  substance  considered.  All 
systems  not  falling  into  a  comparison  on  this  new  basis  have  been 
excluded  from  further  discussion  in  this  report. 


FIGURE  2.  SPECIES  PRESENT  UNDER  EQUILIBRIUM  CONDITIONS  AT  1000  PSIA 
FOR  C/c,  W/c,  AND  TaC/c  EXPOSED  TO  BeF„. 


RHILCQ 

- SSrd^fQtof&m+any, 

AERONUTRONIC  DIVISION 


SECTION  3 

THEORETICAL  RATING  OF  MATERIALS 


3.1  THEORETICAL  RATING  OF  REFRACTORY  MATERIALS  EXPOSED  TO  VARIOUS 
"PURE"  PROPELLANT  SPECIES 

The  theoretical  corrosion  parameters  determined  under  this  con¬ 
tract  for  refractory  materials  exposed  to  various  "pure"  propellant  species 
have  been  correlated  in  terms  of  relative  ratings  of  materials  with  a  given 
species.  The  nominal  species  considered  have  included  AIF3,  BF2,  BF3,  BOF, 
BeF2,  CO,  CO2,  HC1,  HF,  H2,  H2O,  LiF,  N2 ,  and  AI2O3/C.  The  refractories 
C/c,  W/c,  TaC/c  and  HfN/c  were  considered  with  all  14  species,  whereas 
TiC/c,  ZrC/c,  and  were  considered  only  with  the  non-metalllzed 

species,  and  MgO/c  and  Hf02/c,  only  with  C02.  All  comparisons  except 
those  with  AI2O3/C  are  for  a  total  pressure  of  1000  psia.  Graphical 
rating  correlations  for  the  14  esses  are  shown  in  Figures  3  through  16, 
respectively. 

g 

The  primary  rating  parameter  considered  has  been  (volume  of 
refractory  to  saturate  a  unit  weight  of  gas)  as  a  function  of  wall  temper¬ 
ature.  As  previously  stated,  under  assumptions  of  equilibrium  reaction 

g 

with  Lewis  and  Prandtl  numbers  equal  to  unity,  for  a  given  enthalpy 

heat  transfer  coefficient,  should  be  directly  proportional  to  the  linear 
regression  rate  of  the  surface  under  attack.  In  view  of  the  controversial 
aspect  of  calculations  required  to  estimate  over-all  heat  and  mass  transfer 
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FIGURE  3.  THEORETICAL  RATING  OF  REFRACTORY  MATERIALS 
EXPOSED  TO  A1F3  AT  1000  PSIA 
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FIGURE  4.  THEORETICAL  RATING  OF  REFRACTORY  MATERIALS  EXPOSED 
TO  BF2  AT  1000  PSIA 
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FIGURE  5 .  THEORETICAL  RATING  OF  REFRACTORY  MATERIALS  EXPOSED 
TO  BF3  AT  1000  PSIA 
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FIGURE  6 .  THEORETICAL  RATING  OF  REFRACTORY  MATERIALS 
EXPOSED  TO  BOF  AT  1000  PSIA 
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FIGURE  7.  THEORETICAL  RATING  OF  REFRACTORY  MATERIALS 
EXPOSED  TO  BeF ^  AT  1000  PSIA 
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FIGURE  8.  THEORETICAL  RATING  OF  REFRACTORY  MATERIALS 
EXPOSED  TO  CO  AT  1000  PSIA 
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FIGURE  9.  THEORETICAL  RATING  OF  REFRACTORY  MATERIALS  EXPOSED  TO 
C02  AT  1000  PSIA 
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FIGURE  10.  THEORETICAL  RATING  OF  REFRACTORY  MATERIALS 
EXPOSED  TO  HC1  AT  1000  PSIA 
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FIGURE  11.  THEORET ICAL  RATING  OF  REFRACTORY  MATERIALS 
EXPOSED  TO  HF  AT  1000  PSIA 
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(B/p)  ,  CC  WALL  KATEE5L  TO  SATURATE  1  GM  OF 


FIGURE  14.  THEORETICAL  RATING  OF  REFRACTORY  MATERIALS 
EXPOSED  TO  LiF  AT  1000  PSIA 
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FIGURE  15 .  THEORETICAL  RATING  OF  REFRACTORY  MATERIALS 
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EXPOSED  TO  N2  AT  1000  PSIA 
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FIGURE  16.  THEORETICAL  RATING  OF  REFRACTORY  MATERIALS  EXPOSED  TO  Alo0  /c 
AT  THE  MIXED  VAPOR  PRESSURES  OVER  THE  TWO  CONDENSED  PHASES 
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coefficients  for  the  rocket  nozzle  environment,  actual  erosion  rates  them¬ 


selves  have  not  been  predicted  for  the  "pure"  species. 

In  addition  to  the  correlations  of  "isothermal" 


versus  wall 


temperature,  corrosion  heat  effects  also  have  been  considered,  based  upon 
the  "adiabatic  saturation"  model  assuming  an  enthalpy  balance  between 
free  stream  gas,  virgin  wall  material,  and  saturated  mixture  at  the  wall. 
Accordingly,  points  corresponding  to  various  hypothetical  "free  stream" 

g 

temperatures  are  shown  on  the  curves  of  -^5-  versus  wail  temperature  in 
Figures  3  through  16.  The  balance  between  free  stream  and  saturated  wall 
temperatures  does  depend  on  the  assumption  of  unity  Lewis  and  Prandtl 
numbers,  and  hence  lacks  some  of  the  fundamental  significance  of  the 
simple  isothermal  case. 

In  view  of  the  small  variations  of  densities  with  temperature, 
all  values  of  -g-  were  based  upon  room-temperature  densitic  i  for  the  re- 
fractories  in  question.  The  curves  of  versus  wall  temperature  in 
Figures  3  through  16  also  have  been  drawn  as  broken  above  the  respective 
melting  points,  inasmuch  as  the  treatment  assuming  only  chemical  attack 
loses  meaning  where  liquid  run-off  occurs.  (In  rating  materials,  no 
refractory  shall  be  considered  resistant  to  erosion  at  temperatures  above 
its  melting  point.)  Densities  and  melting  points  for  the  nine  refractory 
materials  considered  are  taken  to  be  the  following: 


Refractory 

Density  (929g) 

Melting  Point, 

C/c 

/v  1.90 

.... 

W/c 

19.3 

3650 

TaC/c 

14.5 

4100 

HfN/c 

14.0 

3580 

TiC/c 

4.9 

3410 

ZrC/c 

6.7 

3770 

ZrB2/c 

6.09 

3320 

MgO/c 

3.7 

3070 

Hf02/c 

9.68 

3170 

Relative  ratings  of  refractory  materials  are  discussed  separately  for  the 
fourteen  "pure"  species  considered. 
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AlF-j  -  Figure  3 

Of  Che  four  refractories  studied,  graphite  (C/c)  generally 
appears  the  most  resistant  to  chemical  attack  by  AIF3,  although  TaC/c 
does  give  slightly  lower  values  at  wall  temperatures  from  ^/3800°K 
up  to  its  melting  point  of  4100°K.  Tungsten  (W/c)  also  is  reasonably 
resistant  up  to  its  respective  melting  point.  The  only  material  showing 
serious  theoretical  attack  is  HfN/c,  which  is  oxidized  nearly  quantitatively 
by  the  AIF3 — >A1F  reaction  at  all  temperatures.  Comparisons  on  the  basis 
of  "free  stream"  temperatures  show  almost  the  same  relative  ratings  except 
that  W/c  should  be  cooled  to  below  its  melting  point  of  3650°K  even  at 
bulk  flow  temperatures  well  above  4000°K. 

BF2  -  Figure  4 

All  four  refractories  considered  are  indicated  to  react  nearly 
quantitatively  with  BF2  to  form  borides  at  wall  temperatures  below  about 
2000°K;  only  graphite,  yielding  B^C/c,  exhibits  reasonably  low  values. 
It  is  questionable,  however,  whether  metal  boride  formation  is  kinetically 
likely  at  these  conditions.  Thus,  at  least  in  the  cases  of  TaC/c  and  W/c, 

g 

the  higher  values  at  low  temperatures  might  be  partially  discounted. 

At  high  wall  temperatures,  C/c  and  TaC/c  both  show  excellent 
resistance  to  corrosion  up  to  about  4000°K.  Tungsten  also  is  reasonably 
resistant  up  to  its  melting  point  of  3650°K,  and  should  furthermore  be 
able  to  withstand  "free  stream"  temperatures  of  over  *<000°K  without 
melting. 

Hafnium  nitride  (HfN/c)  reacts  quantitatively  with  BF2  to  form 
HfB2/c  and  Hf-F  species  over  the  entire  temperature  range,  and  thus  cannot 
be  considered  satisfactory  for  use  in  this  atmosphere. 

BF3  -  Figure  5 

Graphite  (C/c)  generally  is  the  most  resistant  to  corrosion  of 
the  four  refractories  considered  with  BF3,  although  all  materials  begin  to 
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show  some  slgnificsnt  attack  at  wall  temperatures  above  3000-3500°K. 

In  the  case  of  tungsten,  in  fact,  the  reaction  to  form  WF^  and  BF  occurs 
nearly  quantitatively  at  its  melting  point  of  3650°K,  Cooling  effects, 
however,  are  significant  with  TaC/c  and  W/c,  so  that  these  two  materials 
show  reasonably  good  corrosion  resistance  at  "free  stream"  temperatures 
up  to  nearly  4000°K. 

Contrary  to  the  case  for  BF2,  metal  boride  formation  with  BF^ 
occurs  only  with  HfN/c,  which  quantitatively  yields  HfBj/c,  Hf-F  species, 
and  correspondingly  high  values  over  the  entire  range  of  temperature. 

BOF  -  Figure  6 

Tungsten  (W/c)  is  the  only  one  of  the  four  refractories  studied 
showing  good  corrosion  resistance  to  BOF.  Above  the  tungsten  melting  point 
of  3650°K,  therefore,  no  material  can  be  rated  as  satisfactory  in  this  en¬ 
vironment  . 

Graphite  (C/c)  and  tantalum  carbide  (TaC/c)  both  react  stolchio- 
metrically  with  BOF  to  form  BF  and  CO  at  wall  temperatures  above  about 
2500°K  and  above  about  3500°K,  respectively.  Reaction  cooling  effects 
are  large,  but  generally  are  still  not  sufficient  to  reduce  corrosion  to 
acceptable  values  at  "free  stream"  temperatures  greater  than  2500°K. 

Hafnium  nitride  (HfN/c)  reacts  quantitatively  with  BOF  over  the  entire 
temperature  range,  yielding  Hf02/c,  HfB2/c  and  Hf-F  species. 

BeF2  -  Figure  7 

Of  the  four  refractories  considered  with  BeFj,  C/c,  W/c  and 
TaC/c  all  show  generally  good  corrosion  resistance  at  wall  temperatures 
up  to  their  melting  points,  or,  in  the  case  of  graphite,  up  to  '''4000°K. 

The  corrosion  reactions  are  mildly  endothermic  in  all  three  cases,  but 
are  not  sufficiently  so  to  provide  cooling  protection  at  "free  stream" 
temperatures  much  above  4000°K.  Except  for  the  case  of  graphite 
volatilisation  at  high  temperatures,  oxidation  by  virtue  of  the 
BeF2  — *BeF  reaction  generally  is  responsible  for  what  attack  does  occur. 
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Hafnium  nitride  (HfN/c)  la  the  only  material  atudied  which  ex¬ 
hibits  appreciable  theoretical  corroaion  by  Be?2  over  the  entire  tern- 
perature  range.  The  values  in  the  vicinity  of  2500°K  wall  temperature 
are  moderately  low,  however,  and  owing  to  the  significant  cooling  effect, 
do  apply  for  "free  stream"  temperatures  up  to  **  3500°K.  Nevertheless, 

HfN/c  must  be  rated  significantly  poorer  than  either  C/c,  W/c  or  TaC/c. 

CO  -  Figure  8 

Of  the  seven  refractory  materials  considered  with  CO,  only  HfN/c 
shows  severe  corrosion  over  the  entire  temperature  range.  The  corrosion 
products  for  HfN/c  and  CO  are  Hf02/c  plus  C/c  below  2550°K,  and  Hf02/c 
plus  HfC/c  above  this  temperature.  The  refractories  TIC/c  and  ZrC/c  are 
oxidized  stoichiometrically  by  CO,  but  only  at  wall  temperatures  of 
2000°K  and  below.  At  2500°K,  the  oxidation  becomes  nil.  Tungsten  also 
shows  corrosion  by  virtue  of  the  formation  of  WC/c  at  low  wall  temperatures, 
but  this  reaction  may  be  kinetically  unlikely  at  these  conditions. 

Other  than  the  preceding,  all  materials  considered,  including 
TaC/c,  ZrBj/c ,  and  C/c,  show  excellent  resistance  to  attack  by  CO  at 
temperatures  up  to  their  melting  points,  or,  in  the  case  of  C/c,  until 
simple  volatilization  becomes  important. 

CO2  -  Figure  9 

Of  the  nine  refractory  materials  considered  with  CO2,  only  the 
oxides  MgO/c  and  Hf02/c  show  some  satisfactory  resistance  to  chemical 
corrosion.  Use  of  these  two  refractories,  however,  must  take  into  account 
their  questionable  physical  properties  and  their  low  melting  points  (3070°K 
and  3170°K,  respectively). 

The  remaining  seven  refractories  are  all  attacked  nearly  stoichio¬ 
metrically  by  CO2,  forming  CO  in  the  case  of  C/c,  gaseous  W3O9  in  the  case 
of  W/c,  and  condensed  metal  oxides  with  the  others.  Over  the  range  of 
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wall  temperature*  considered,  tungsten  yields  slightly  lower  values 
than  the  other  materials,  although  if  "free  stream"  temperature  is  taken 
as  the  variable,  graphite  is  more  resistant  below  2500°K,  owing  to  a 
pronounced  cooling  effect. 

HC1  -  Figure  10 

Tungsten  and  graphite  both  exhibit  excellent  resistance  to 
corrosion  by  HC1  at  wall  temperatures  up  to  3650°K,  the  melting  point  of 
tungsten.  Tantalum  carbide  (TaC/c)  becomes  increasingly  resistant  to 
attack  at  the  higher  wall  temperatures,  and,  between  3650°K  and  its  own 
melting  point  of  4100°K,  appears  to  be  the  best  choice.  On  the  basis  of 
"free  stream"  temperature,  however,  TaC/c  loses  its  melting  point  advantage 
over  W/c,  owing  to  an  exothermic  heat  of  reaction. 

The  remaining  four  refractories  considered  with  HC1  (HfN/c,  TIC/c, 
ZrC/c  and  ZrB2/c)  all  show  moderate  to  severe  theoretical  corrosion  over 
the  entire  temperature  range,  resulting  from  near- stoichiometric  reactions 
forming  gaseous  metal  chlorides  and  hydrogen. 

HF  -  Figure  11 

Only  graphite,  of  the  seven  refractories  considered,  exhibits 
good  resistance  to  corrosion  by  HF,  and  only  at  wall  temperatures  below 
3500°K  ("free  stream"  temperatures  below  4000°K) .  Tungsten  and  tantalum 
carbide  are  attacked  nearly  stoichiometrically  by  HF  at  low  temperatures, 
but  are  somewhat  more  resistant  at  higher  temperatures,  so  that  they  may 
be  rated  roughly  equivalent  to  graphite  at  temperatures  in  the  vicinities 
of  their  melting  points. 

The  four  remaining  refractories  HfN/c,  TIC/c,  ZrC/c  and  Zrl^/c 
undergo  nearly  quantitative  reactions  with  HF  throughout  most  of  the  tem¬ 
perature  range,  yielding  metal  fluorides  and  hydrogen. 
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H2  -  Figure  12 

The  only  refractory  material  of  the  aeven  conaldered  that  la 
aeverely  attacked  by  H2  under  equilibrium  condition!  la  graphite  (C/c). 

At  temperature!  below  2000°K  the  primary  reaction  product  la  CH^;  above 
2500°K,  the  products  are  C2H2  and  C2H.  The  hydrogenolysls  of  graphite  la 
recognized  to  be  somewhat  kinetlcally  limited  in  practice. 

The  other  six  materials  studied,  W/c,  TaC/c,  HfN/c,  TIC/c, 

ZrC/c  and  ZrB2/c,  all  show  excellent  resistance  to  corrosion  by  H2  at 
temperatures  up  to  their  respective  melting  points.  Condensed  metal 
hydrides,  however,  particularly  HfH2/c,  ZrH2/c  and  T1H2/c,  were  not  con¬ 
sidered  in  the  calculations  from  which  the  curves  In  Figure  12  were 
derived. 

H20  -  Figure  13 

All  seven  refractories  considered  show  severe  theoretical 
corrosion  by  HjO.  The  most  nearly  resistant  material  over  the  entire 
temperature  range  is  tungsten  (W/c) ,  which  is  also  the  only  refractory 
not  oxidized  quantitatively  at  these  conditions.  Nevertheless,  values 
for  tungsten  are  high.  The  remaining  six  materials  are  converted  to  CO, 
in  the  case  of  C/c,  and  condensed  metal  oxides  in  the  cases  of  TaC/c, 

HfN/c,  TiC/c,  ZrC/c  and  ZrB2/c. 

L1F  -  Figure  14 

At  temperatures  much  below  3000°K,  systems  consisting  only  of  a 
refractory  material  and  pure  LiF  are  completely  condensed  at  1000  psla. 
Restricting  attention  to  the  higher  temperature  regime,  TaC/c  and  W/c 
show  excellent  resistance  to  corrosion  by  LiF  up  to  their  respective 
melting  points  (4100°K  and  3650°K) .  Graphite  also  is  quite  resistant  up 
to  the  temperature  (4000-4500°K)  at  which  vaporization  becomes  important. 

Of  the  four  materials  considered,  only  HfN/c  exhibits  appreciable  theoretical 
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corrosion  at  elevated  wall  temperatures.  Even  then,  owing  to  the  pro* 
nounced  cooling  effect,  HfN/c  may  be  moderately  realatant  to  attack  If 
"free  atream"  temperature  la  taken  aa  the  baala. 


N2  -  Figure  15 

Of  the  aeven  refractorlea  conaldered  with  a  N2  environment, 

TIC/c,  ZrB2/c  and  poaalbly  TaC/c  and  ZrC/c  exhibit  aevere  theoretlcel 
corroalon  at  temperaturea  in  the  range  1500-2500°K,  owing  to  the  stolchlo- 
metric  formation  of  metal  nitridea.  Aa  in  the  caaea  of  boride  and  carbide 
formation  at  low  temperatures,  the  formation  of  nitrides  at  auch  conditions 
may,  however,  be  kinetically  unlikely. 

At  higher  wall  temperaturea,  nearly  all  the  materials  conaldered 
were  quite  resistant  to  corrosion  by  N2  up  to  their  melting  points;  the 

g 

very  lowest  values  were  found  for  W/c  and  for  HfN/c .  Only  In  the  case 
of  graphite  did  theoretical  attack  begin  to  become  appreciable  above  3500°K, 
owing  to  the  formation  of  gaseous  CN. 

AI2O3/C  -  Figure  16 

Inasmuch  as  most  systems  of  a  refractory  material  plus  AI2O3/C 
would  be  completely  condensed  at  1000  psla  except  at  very  high  temperatures, 
a  comparison  of  equilibrium  corrosion  behavior  for  a  fixed  total  pressure 
was  not  considered  meaningful.  Instead,  calculations  for  AI2O3/C  and  the 
four  refractories  C/c,  W/c,  TaC/c,  and  HfN/c  were  carried  out  for  variable 
total  pressures,  such  that  the  two  condensed  phases,  AI2O3/C  and  the  re¬ 
fractory,  would  both  be  present  at  equilibrium.  From  phase-rule  con¬ 
siderations,  only  one  such  mixed  "vapor  pressure"  and  one  gas  composition 

could  exist  for  any  given  temperature. 

B 

Comparison  of  curves  in  Figure  16  Indicates  that  W/c  is  the 
most  resistant  to  attack  by  AI2O3/C  at  these  conditions,  yielding  only 
moderate  theoretical  corrosion  up  to  its  melting  point  of  3650°K.  In 
contrast,  C/c,  TaC/c  and  HfN/c  all  attacked  almost  quantitatively  by 
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AljOj/c,  yielding  CO  in  Che  ceae  of  graphite  end  TA/c  and  HfC^/c  in  the 
caae  of  the  metallic  refractoriea .  Vapor  preaaurea  over  the  two  phaaea 
alao  were  loweat  for  W/c  plua  AI2O3/C;  a  total  preaaure  of  1000  pala  would 
not  be  reached  except  at  a  temperature  above  5000°K,  i.e.,  much  higher 
than  the  melting  point  of  tungaten.  Total  preaaurea  of  1000  pala  would 
be  reached  inatead  at  about  2900°K  for  C/c,  3800°K  for  TaC/c,  and  3500°K 
for  HfN/c,  At  temperaturea  higher  than  the  preceding,  "boiling"  would 
take  place  if  the  two  condenaed  phaaea  were  in  contact  at  1000  pala  total 
preaaure . 

Conclusion 

Relative  theoretical  ratings  for  refractory  materials  exposed 
to  the  fourteen  propellant  species  studied  are  sumnarlzed  in  Table  III. 

It  is  seen  that  with  only  a  few  exceptions,  the  beat  policy  from  the 
standpoint  of  corrosion  apparently  fits  the  empirical  rule  of  choosing 
tungsten  for  wall  temperatures  up  to  its  melting  point  (3650°K) ,  and 
then  choosing  graphite  for  higher  wall  temperatures. 
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RELATIVE  THEORETICAL  RATINGS  OF  REFRACTORY  MATERIALS 
EXPOSED  TO  14  DIFFERENT  PROPELLANT  SPECIES 


C/c 

W/c 

TaC/c 

HfN/c  : 

TiC/c 

ZrC/c  ZrB2/c 

MgO/c 

Hf02/c 

A1F3 

® 

X 

X 

0 

- 

- 

- 

- 

bf2 

® 

? 

? 

0 

- 

- 

- 

- 

- 

BF3 

X 

X 

0 

- 

- 

- 

- 

- 

bof 

0 

0 

0 

0 

- 

- 

- 

- 

- 

BeF2 

X 

X 

® 

0 

- 

- 

- 

- 

- 

CO 

X 
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® 

0 

0 

0 

X 

- 

- 

co2 

0 

•o') 

v  ^ 

0 

0 

0 

0 

0 

X 

® 

HC1 

X 

® 

X 

0 

0 

0 

0 

“ 

- 

HF 

0 

0 

0 

0 

0 

0 

- 

- 

H2 

0 

X 

X 

1 

i  x 

X 

X 

- 

- 

h2o 

0 

o » 

0 

0 

0 

0 

0 

- 

- 

L1F 

X 

X 

® 

0 

- 

- 

- 

- 

- 

N2 

X 

® 

? 

X 

? 

? 

? 

- 

- 

A1203/c 

0  ®  0  0 

KEY:  X  "acceptable" 

?  questionable 

0  not  "acceptable" 

-  not  studied 

(Best  choices  are  circled  for 

each 

species.) 
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3.2  THEORETICAL  RATING  OF  REFRACTORY  MATERIALS  EXPOSED  TO  REAL  PROPELLANT 
SYSTEMS 

In  addition  to  the  atudlea  of  theoretical  corroaion  of  refractory 
materiala  by  "pure"  combustion  apeciea,  a  few  ainilar  calculationa  were 
made  for  materiala  expoaed  to  real  propellant  mixturea.  The  comparisons 
to  be  diacuaaed  in  the  following  paragrapha  include  auch  ayatema  ap  F2/H2, 
aluminized  "hot  aolid",  and  "NF"-"BOF"  type  combinationa. 

F2/H2  -  Figure  17 

Theoretical  aaturation  parameters  (—*—)  for  four  refractoriea 
(TIC/c,  TiN/c,  ZrC/c  and  Zr02/c)  expoaed  to  a  F2/H2  ayatem  at  500  psia  are 
compared  in  Figure  17.  The  principal  species  present  in  the  "free  stream" 
combustion  mixture  at  the  throat  are  HF,  H2  and  H.  All  the  materials  con¬ 
sidered  showed  appreciable  theoretical  corrosion  through  the  temperature 
region  of  interest.  Considering  that  Zr02/c  would  also  undergo  melting 
ablation  at  wall  temperatures  above  2950°K,  perhaps  the  moat  resistant  of 
the  four  would  be  TIC/c.  However,  graphite  and  (possibly)  tungsten, 
neither  of  which  were  considered  with  this  particular  mixture,  would  both 
probably  be  rated  more  satisfactory  than  TiC/c. 

"Hot  Solid"  -  Figures  18  and  19 

Theoretical  saturation  values  are  plotted  in  Figure  18  for  C/c, 
W/c  and  TaC/c  refractories  exposed  to  the  combustion  products  of  an  alum¬ 
inized  solid  propellant  containing  a  perchlorate  oxidizer.  The  calculations 
are  for  a  throat  pressure  of  405.8  psia,  corresponding  to  a  chamber  pressure 
of  700  psia.  The  major  species  present  in  the  "free  stream"  at  the  throat 
are  CO,  H2O,  H2,  H,  N2,  HC1  and  AI2O3/C.  Under  these  conditions,  tungsten 
(W/c)  is  clearly  the  moat  resistant  to  attack  at  wall  temperatures  up  to 
its  melting  point  of  3650°K.  Since  tungsten,  however,  would  undergo 
melting  ablation  at  the  wall  temperature  corresponding  to  the  3826°K 
throat  stagnation  temperature,  none  of  the  three  materials  can  be  rated 
completely  satisfactory  for  this  service. 
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FIGURE  17.  THEORETICAL  RATING  OF  REFRACTORY  MATERIALS 
EXPOSED  TO  F2/H2  (0/F  -  5.50)  AT  500  PSIA 
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FIGURE  18.  THEORETICAL  RATING  OF  REFRACTORY  MATERIALS  EXPOSED  TO 

"HOT  SOLID"  PROPELLANT*  AT  405.8  PSIA  THROAT  PRESSURE** 
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FIGURE  19.  THEORETICAL  RATING  OF  REFRACTORY  MATERIALS  EXPOSED  TO 
"HOT  SOLID"  PROPELLANT  (MINUS  Al  0  )*  AT  405.8  PSIA 
THROAT  PRESSURE**  1  ' 


FIGURE  20.  THEORETICAL  RATING  OF  REFRACTORY  MATERIALS  EXPOSED  TO 
"NF"  PROPELLANT*  AT  287  PSIA  THROAT  PRESSURE** 
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FIGURE  21.  THEORETICAL  RATING  OF  REFRACTORY  MATERIALS  EXPOSED  TO 

SIMULATED  "BOF"  PROPELLANT*  AT  229.63  PSIA  THROAT  PRESSURE** 
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In  Che  caee  of  e  solid  propellent  yielding  condensed  AljO^ 
particles,  theoretical  corrosion  calculations  can  be  made  either  for  the 
assumption  that  the  particles  react  with  the  wall  in  the  same  proportion 
as  their  bulk  gas  concentration,  as  in  the  preceding  example,  or  for  the 
alternative  assumption  that  the  particles  remain  in  the  turbulent  core 
and  never  penetrate  the  boundary  layer.  In  this  second  case,  which  is 
shown  in  Figure  19  for  the  three  refractories  C/c,  W/c  and  TaC/c,  the 
Al^j/c  present  in  the  mixture  at  throat  stagnation  conditions  has  been 
"removed"  from  the  system  prior  to  equilibration  with  the  wall  material. 

It  is  seen  that  this  assumption  yields  only  minor  changes  in  absolute 
values  versus  wall  temperatures,  and  no  changes  at  all  in  the  relative 
ratings.  In  other  words,  the  corrosivity  of  the  gaseous  phase  is  of 
comparable  magnitude  to  that  of  the  Al^O^/c. 

Calculations  of  theoretical  corrosion  for  an  approximate  simu¬ 
lation  of  the  "hot  solid"  propellant  at  a  slightly  different  pressure 
(not  shown)  indicated  only  small  differences  from  the  results  given  in 
Figure  18  and  19. 

"NF"-"BOF"  Types  -  Figures  20  and  21 

Equilibrium  corrosion  parameters  for  C/c,  W/c  and  TaC/c  re¬ 
fractories  exposed  to  an  "NF"  solid  propellant  at  287  psla  throat  pressure 
(500  psla  chamber  pressure)  are  shown  in  Figure  20.  The  major  combustion 
products 'of  this  composition  at  the  throat  conditions  are  CO2,  CO,  H20,  H2, 
H,  N2,  HF  and  BOF .  Tungsten  once  again  is  seen  to  be  the  most  resistant 
to  attack,  but  owing  to  melting  ablation,  still  is  not  capable  of  with¬ 
standing  the  3911°K  throat  stagnation  temperature.  The  relative  ranking 
of  TaC/c  and  C/c  varies,  depending  on  whether  TaN/c  formation  is  assumed 
to  be  kinetlcally  likely  with  the  former,  but  in  any  case,  attack  should 
be  appreciable  with  either  of  these  two  materials. 

In  Figure  21  the  three  materials  C/c,  W/c  and  TaC/c  again  are 
compared,  this  time  for  exposure  to  an  approximate  simulation  of  a  "BOP" 
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propellant  at  229.63  pala  throat  preaaure  (400  pala  chamber  preaaura). 

The  major  epeclea  here  are  CO,  H2O,  H2,  H,  N2,  HF  and  BOF.  Inaemuch  aa 
the  chemical  coopoaition  la  not  too  different  from  that  of  the  "HF" 
propellant,  the  relative  ratlnga  are  much  the  same,  with  W/c  ahowing  ex¬ 
cellent  corroalon  reaiatance  at  wall  temperaturea  up  to  lta  melting  point, 
while  C/c  and  TaC/c  are  attacked  aubatantlally  throughout  the  temperature 
regime.  (The  case  of  TaN/c  formation  auppreased  was  not  considered  for 
TaC/c  In  this  comparison.) 
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SECTION  4 

PREDICTION  QF  ACTUAL  EROSION  RATES 


In  the  Fifth  Quarterly  Report  under  thi8  contract,  a  pre¬ 
liminary  study  was  presented  which  attempted  to  assess  the  validity  of 
he"simplif led  Denison"  ablation  model^assuming  equilibrium  at  the  wall 
with  Lewis  and  Prandtl  numbers  taken  to  be  unity.  A  three-way  comparison 
was  made  between  predicted  erosion  rates,  rates  calculated  by  Jones  and 
Delaney®  from  kinetic  theory,  and  rates  determined  experimentally  by 

q 

NASA-Lewis.  The  systems  considered  were  graphitic  nozzle  Inserts  ex¬ 
posed  to  two  state-of-the-art  solid  propellants  and  to  two  advanced 
solid  formulations,  with  nominal  chamber  pressures  of  1000  psia  for 
the  former  and  500  psia  for  the  latter. 

The  erosion  rates  predicted  in  this  study  were  estimated  from 
equation  (2)  or  equation  (3)  for  the  unsteady-state  simplified  Denison 
model  with  wall  temperature  assumed  to  be  known  as  a  function  of  time. 

Owing  to  the  relatively  short  duration  of  the  firings  being  considered, 

the  steady-state  adiabatic  saturation  model  (equation  (1)  )  was  not 

^c 

applicable.  The  enthalpy  heat  transfer  coefficient  (gjj-)  was  estimated 
using  the  Bartz ^equation  for  rocket  nozzle  heat  transfer,  even  though 
this  correlation  was  recognized  to  be  only  an  approximation.  For  those 
propellants  yielding  condensed  AI2O3/C  as  a  combustion  product,  theoretical 
corrosion  calculations  were  performed  both  for  the  assumption  that  the 
AI2O3/C  particles  diffuse  to  the  wall  in  proportion  to  their  concentration, 
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and  for  the  alternative  assumption  that  no  AI2O3/C  ever  reaches  the  wall. 
(Strictly  speaking,  the  Denison  approachable  limited  to  all -gas  systems, 
and  thus  should  be  applicable  only  to  the  latter  case.) 

The  experimental  data  considered  were  determined  at  MASA-Lewls^ 
for  three  types  of  graphite  nozzles  (ATJ,  Speer  3499,  and  ZT)  exposed  to 
the  combustion  products  of  Arcite  368  and  Arcite  373,  two  state-of-the-art 
solid  propellants.  The  former  of  tnese  is  a  4700°F  nonmetallized  composi¬ 
tion  containing  the  elements  CH0NC1,  and  the  latter,  a  5500°F  aluminized 
formulation  including  the  elements  CH0NC1A1.  Inasmuch  as  the  unpublished 
results^for  test  firings  with  Arcite  373  showed  erosion  completely  ob¬ 
scured  by  deposition  of  a  layer  of  AI2O3/C,  only  the  data  for  Arcite  368 
were  considered  for  comparison. 

Predictions  of  theoretical  corrosion,  using  the  simplified 
Denison  model  assuming  equilibrium  reaction,  averaged  almost  a  full  order 
of  magnitude  higher  (~35  mils /second)  than  the  maximum  rates  calculated 
(~5  mils/second)  for  the  NASA  tests  with  Arcite  368  and  graphite  nozzles. 
The  discrepancy  was  explained  by  examination  of  surface  temperatures,  which, 
as  estimated  by  a  somewhat  doubtful  extrapolation  of  thermocouple  readings 
for  the  ATJ  firing,  never  exceeded  1800°K  during  the  30-second  duration 
of  the  test.  At  temperatures  this  low,  surface  reaction  rates  might  well 
be  nil,  rather  than  infinitely  rapid  as  was  assumed  in  the  equilibrium 
model.  From  chamber  pressure  traces  taken  during  the  NASA  tests,  it 
appears  that  most  of  the  measured  erosion  of  graphitic  nozzles  was  purely 
mechanical  in  nature,  and  occurred  during  the  first  five  seconds  of  firing. 

Predictions  based  upon  the  simplified  Denison  approacu  also  were 
compared  to  erosion  rates  calculated  by  Jones  and  Delaney**  for  graphite 
nozzles  exposed  to  various  propellant  systems.  In  most  cases,  the  cal¬ 
culations  of  Jones  and  Delaney  were  based  upon  finite  reaction  rates, 

-  n 

rather  than  attainment  of  complete  chemical  equilibrium.  Independent 
correlations  for  overall  heat  transfer  and  for  overall  mass  transfer 
moreover  were  employed,  so  that  the  assumption  of  unity  Lewis  numbers  was 
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not  required.  For  these  reasons,  Che  predictions  of  Jones  and  Delaney 
would  be  expected  to  differ  from  those  based  upon  the  simple  model, 
and  should  probably  be  better. 

In  the  case  of  the  Arclte  368  propellant,  and  for  comparable 
conditions  of  temperature,  pressure,  and  geometry,  the"slmpllf led  Denison" 
predictions  assuming  equilibrium  corrosion  averaged  an  order  of  magnitude 
higher  than  the  erosion  rates  calculated  by  Jones  and  Delaney  for  kinetic- 
ally-limlted  reaction.  The  two  methods,  however,  were  In  rough  agreement 
(factor  of  only  1.5 -1.7)  for  the  special  case  In  which  Jones  and  Delaney 
also  assumed  complete  equilibrium  at  the  wall.  Neither  method  predicted, 
or  could  predict,  the  mechanical  erosion  effect  previously  described 
,  as  having  occurred  in  the  NASA  experiments  with  this  system. 

No  comparison  with  meaningful  test  data  was  possible  for 
Arclte  373  propellant,  nor  for  the  two  advanced  formulations,  "NP"*  and 
"NF"**,  respectively.  Comparisons  were  made  between"simplif led  Denison" 
predictions  and  the  more  rigorous  calculations  of  Jones  and  Delaney  for 
these  systems ,  however,  and  showed  graphite  erosion  rates  estimated  by 
the  simple  equilibrium  model  high  by  a  factor  of  at  least  3  for  Arclte  368 
and  for  "NP",  and  high  by  a  factor  of  about  8  for  the  "NF"  case.  Inasmuch 
as  kinetic  limitations  probably  become  unimportant  at  the  high  (/^3000°K) 
wall  temperatures  involved,  the  discrepancies  between  the  two  methods  must 
arise  either  from  (1)  failure  of  the  Denison  approach  to  allow  for  Lewis 
number  effects,  or  from  (2)  omission  of  important  species  or  reactions  in 
the  calculations  by  Jones  and  Delaney.  (The  BOF-graphite  reaction,  for 
example,  was  neglected  in  the  "NF"  propellant  case.) 

In  conclusion,  the  simplified  Denison-model  approach  assuming 
chemical  equilibrium  at  the  refractory  wall  is  not  valid  for  predicting 


*  "NP"  *  6600°F  solid  propellant,  elements  CHONCILiAl 
**"NF"  »  6700°F  solid  propellant,  elements  CHONFB 
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transient  erosion  rates  with  most  state-of-the-art  solid  propellant 
svstems,  owing  to  the  low  temperatures  involved.  Por  advanced  solid 
propellants  yielding  higher  flame  temperatures,  the  simple  method  does 
show  some  promise,  at  least  for  making  qualitative  comparisons  for 
different  refractories.  Greater  quantitative  validity  for  high-tempera* 
ture  systems  might  be  achieved  if  the  simplified  Denison  model  could  be 
modified,  particularly  to  account  for  the  effects  of  non-unity  Lewis 
numbers . 


I 

I 

1 
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SECTION  5 

PROPELLANT  PERFORMANCE  STUDIES 


Theoretical  performance  calculations  for  certain  unclassified 
and  classified  binary  liquid  propellant  combinations  were  carried  out  as 
part  of  this  contract.  The  "pressure  profile"  option  of  the  Aeronutronic 
Rocket  Propellant  Performance  Program  was  utilized  in  determining  expansion 
performance  parameters  at  integer  area  ratios  from  1  to  50,  and,  in  turn, 
ideal  performance  values  at  given  epsilons  for  ambient  pressures  ranging 
from  sea  level  down  to  vacuum.  Thus,  propellant  systems  could  be  compared 
for  specific  applications,  rather  than  on  some  arbitrary  basis  such  as 
optimum  expansion  to  one  atmosphere. 

Calculations  were  run  for  two  chamber  pressures:  1000  psia  and 
300  psia.  In  general,  five  mixture  ratios  at  each  pressure  were  considered 
for  each  propellant  combination  studied.  The  latest  JANAF3-^ thermodynamic 
data  available  at  the  time  were  used,  with  a  large  number  of  possible 
combustion  species  Included.  Both  the  frozen  and  the  equilibrium  expansion 
cases  were  considered. 

The  results  for  eleven  unclassified  systems  (CIF3/N2H4,  F2/H2, 

H202/B5H9,  n2f4/b5h9,  n2h4/b5h9,  n2o4/n2h4,  of2/b2h6,  of2/b5h9,  of2/n2h4, 

O2/H2,  and  O2/RP-I)  have  been  reported  in  Volume  II  of  the  Second  Quarterly 
Report6.  Similarly,  the  data  for  ten  classified  systems  have  been  presented 
in  Volume  III^  of  that  report.  Inasmuch  as  the  method,  computer  program, 
and  results  were  discussed  fully  at  that  time,  the  reader  is  referred  to 
these  two  volumes  for  further  information. 
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APPENDIX  A 

THERMODYNAMIC  DATA  FOR  C2H,  THREE  W 
SPECIES,  AND  FIVE  Zr  SPECIES 
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37.587 

44.587 

1600.00 

133.259 

163.476 

37.602 

48.347 

1700.00 

135.105 

165.756 

37.614 

52.107 

1800.00 

136.868 

167.906 

37.624 

5S.869 

1900.00 

138.555 

169.941 

37.633 

59.652 

2000.00 

140.173 

171.871 

37.640 

63.396 

2100.00 

141.727 

173.708 

37.646 

67.160 

2200.00 

143.221 

175.459 

37.6S2 

70.925 

2300.00 

144.659 

177. 133 

37.656 

74.690 

2400.00 

146.046 

178.736 

37.661 

78.456 

2500.00 

167.384 

100.273 

37.664 

82.222 

2600.00 

148.678 

181.751 

37.668 

85.989 

2700.00 

149.929 

183.172 

37.671 

89.756 

2800.00 

151.  141 

184.542 

37.673 

93.523 

2900.00 

152.316 

185.864 

37.676 

97.291 

3000.00 

153.455 

107. 142 

37.678 

101.050 

3100.00 

154.562 

108.377 

37.680 

104.826 

3200.00 

155.638 

189.573 

37.601 

108.594 

3300.00 

156.684 

190.733 

37.683 

112.562 

3400.00 

157.702 

191.058 

37.684 

116.131 

3500.00 

158.693 

192.950 

37.686 

1 19.099 

3600.00 

159.660 

194.012 

37.687 

123.668 

3700.00 

160.602 

195.044 

37.608 

127.437 

3000.00 

161.522 

196.050 

37.689 

131.205 

3900.00 

162.420 

197.029 

37.690 

134.974 

4000.00 

163.297 

197.983 

37.691 

138.744 

6100.00 

164.154 

198.914 

37.692 

142.513 

4200.00 

164.993 

199.022 

37.693 

146.282 

6300.00 

165.813 

200.709 

37.693 

150.051 

4400.00 

166.616 

201.575 

37.694 

153.821 

4500.00 

167.402 

202.422 

37.695 

157.590 

4600.00 

168. 173 

203.251 

37.695 

161.359 

4700*00 

168. 928 

204.062 

37.696 

165. 129 

6000.00 

169.660 

204.855 

37.696 

168.899 

4900.00 

170.394 

205.632 

57.697 

172.668 

5000.00 

171. 106 

206.394 

37.697 

176.438 

5100.00 

171.806 

207. 141 

37.698 

100.208 

5200.00 

172.492 

207.873 

37.698 

183.977 

5300.00 

173. 167 

208.591 

37.698 

187.747 

5400.00 

173.829 

209.295 

37.699 

191.517 

5500.00 

174.480 

209.987 

37.699 

195.287 

5600.00 

175.120 

210.666 

37.699 

199.057 

5700.00 

175.750 

211.354 

57.700 

202.027 

5800.00 

176.369 

211.989 

37.700 

206.597 

5900.00 

1  76.978 

212.634 

37.700 

210.367 

6000.00 

177.578 

213.267 

57.700 

214.137 
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THBueonuMic  ruicnoNS  or  t»(») 
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T 

(°K> 

)/T 

(cal/aol*  °l) 

S 

(cal/aol*  °K) 

CP 

(ctl/aPli  °K> 

HT 

“2*8 

(kcal/aoU) 

2»8.  16 

5M.9C0 

5H.900 

M.090 

0 

300.00 

56.900 

68.950 

8.097 

.016 

••CO. 00 

69.222 

6 1 . 325 

8.400 

.841 

500.00 

57. Hid 

61.220 

8.572 

1.691 

600.00 

60.617 

64. 792 

8.676 

2.553 

700.00 

61.241 

66. 115 

8.742 

3.425 

600.00 

61.929 

67.105 

8.7  87 

4.301 

700.00 

02.685 

68. 142 

8.818 

5.  181 

1000.00 

6  1.206 

69.271 

8.841 

6.064 

1  100.00 

63. 798 

70.  1  16 

8.858 

6.949 

1200.30 

64.167 

70.887 

8.871 

7.M36 

1100.00 

64.887 

7 1 . 598 

8.882 

8.724 

1400.00 

65.390 

72.256 

8.890 

9.612 

1500.00 

65.867 

72.870 

8.897 

10.502 

1600.00 

66.325 

71.444 

8.902 

11.391 

1700.00 

66.757 

71.984 

8.907 

12.282 

1600.00 

67.  1  75 

74.493 

8.910 

13.173 

1900. OC 

67.571 

74.975 

8.914 

14.064 

2000.00 

67.955 

75.412 

8.916 

14.956 

2100.00 

68.321 

75.867 

8.919 

15.847 

2200.00 

68.674 

76.282 

8.921 

16.739 

2100.00 

69.013 

76.679 

8.923 

17.631 

2400.00 

69.341 

77.059 

8.924 

18.524 

2500.00 

69.657 

77.423 

tt.  926 

IV. 416 

2600.00 

69.962 

- T-T-i-7  7 1 

8.927 

20.309 

2700.00 

70.258 

78.  1  10 

8.928 

21.202 

2800.00 

70.544 

78.416 

8.929 

22.095 

2900.00 

70.821 

78. 748 

8.930 

22.988 

1000.00 

71.091 

79.051 

8.931 

23.881 

1100.00 

71.352 

79.344 

8.932 

24.774 

3200.00 

71.606 

79.627 

8.932 

25.667 

3100.00 

71.854 

79.902 

0.913 

26.560 

3400.00 

72.094 

80. 169 

8.911 

27.453 

3500. OC 

72.329 

80.428 

8.934 

28.347 

3600.00 

72.567 

80.680 

8.934 

29.240 

3700.00 

72.780 

80.924 

8.935 

30.134 

1800.00 

72.998 

81.  163 

8.936 

31.027 

3900.00 

71.210 

81.395 

8.936 

31.921 

4000.00 

71.417 

81.621 

8.936 

32.814 

4  100.00 

71.620 

81.842 

8.936 

33.708 

4200.00 

73.818 

82.057 

8.916 

34.602 

4300.00 

74.011 

82.267 

8.937 

35.495 

4400.00 

74.202 

82.473 

8.937 

36.389 

4500.00 

74.189 

82.674 

6.937 

37.283 

4600.00 

74.571 

82.870 

8.93  7 

38.176 

4700.00 

74.749 

81.062 

8.938 

39.070 

4800.00 

74.925 

83.250 

8.938 

39.964 

4900. 00 

75.096 

83.415 

8.938 

40.858 

5000.00 

75.265 

83.615 

8.938 

41.751 

5100.00 

75.410 

83.792 

8.938 

42.645 

5200.00 

75.593 

81.966 

8.939 

43.539 

5100.00 

75.752 

84. 116 

8.939 

44.433 

5400.00 

75.909 

84.101 

8.939 

45.327 

5500. OC 

76.061 

84.467 

8.939 

46.221 

5600.00 

76.215 

84.628 

8.939 

47.115 

5700.00 

76.164 

84.786 

8.939 

48.009 

5800.00 

76.510 

84.942 

8.939 

48.902 

5900.00 

76.655 

85.095 

8.919 

49.796 

6000.00 

76.797 

85.246 

8.939 

50.690 
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nuMonumc  rnmiM  or  t»«}/c 

A»°jm  >  -74.3  (CU/KU 

^  s  "m 

(c«l/Mli  JO  °K) _ (ctl/melt  TO  (kc«l/«oU) 


t 

(°K) 


296.13 
100.00 
400.00 
100.00 
tco.oc 
700.00 
800.00 
900. 0C 
IOCO.OO 
1100.00 
1200.00 
1100.00 
uco.oo 

1500.00 
1600.00 
1700.00 
1600.00 
1900.00 
2000. OC 
2  ICO. 00 
2200.00 
2100.00 
2400. OC 
2100.00 
2600.00 
2700.00 
28CO.CC 
2900.00 
3000. OC 
1100.00 
1200.00 
1100.00 
3400.00 
1100.00 
1600.00 
1700. OC 
38CO.OC 
1900.00 
4000. OC 
4100.00 
42CC.0C 
4100.00 
4400.00 
4100.00 
460C.CC 
4700.00 
4800.00 
4900. OC 
IOCO.OO 
1100.00 
1200. OC 
1300. OC 
14CC.0C 
1100.00 
S6CO.CO 
17C0.0C 
38C0.CC 
1900.00 
6000. OC 


8.190 

8.190 

9.111 

10.101 

11.664 

13.081 

14.491 

11.864 

17.181 

18.411 
19.661 
20.622 

21.932 
22.9V6 
24.016 
24.996 
21.939 
26.847 
27.721 
28.169 
29.386 
30.178 
10.941 
31.690 

32.412 
11.111 
11.798 
14.464 

31.111 
11. 746 
16.361 
36.966 
37.660 
18.149 
19.019 
39.670 
4C.103 
40.920 

41.121 
42. 107 
42.678 
43.236 
43.782 
44.314 
44.831 
41.341 
41.844 
46.332 
46.811 
47.280 
4  7.  74C 
48.191 
46.633 
49.068 
49.494 

49.913 
10.321 
10.729 
11.127 


8.190 

8.471 

12.913 
16.766 
20.101 
21.032 

21.640 
27.987 

30.121 
32.078 
33.886 
31.167 
17.  140 
18.617 
40.012 
41.311 
42.189 
41.787 

44.932 
46.029 
47.081 
48.096 
49.077 
10.022 
10.917 
1I.H23 
12.681 
13.118 
14.331 

11.122 

11.895 

16.641 
61.431 
62.116 
62.778 
63.422 
64.048 
64.619 
61.214 
65.614 
66.400 

66.913 
67.494 
68.022 

68.138 
69.044 

69.138 
70.023 
70.498 
70.963 
71.419 
71.867 
72.306 
72.737 

73.161 
71.177 
71.981 
74.367 
74.782 


13.112 
13.198 
16.294 
17.610 
18.704 
19.306 
19.730 
20.104 
20.401 
20.662 
20.897 
21.114 
21.319 
21.114 
21.702 
21.881 
22.061 
22.217 
22.409 
22.178 
22.746 
22.912 
23.077 
21.240 
21.403 

21.161 
21.726 
21.887 
24.047 
24.206 
24.366 

24.121 

23.100 

23.100 

23.100 
21. IOC 
21. IOC 

23.100 

23.100 

21.100 

21.100 

21.100 

21.100 

21.100 

21.100 

21.100 

21.100 

23.100 

21.100 

21.100 

21.100 

21.100 

21.100 

21.100 

21.100 

23.100 

21.100 

21.100 

21.100 


0 

.024 

1.120 

3.211 

1.062 

6.964 

8.918 

10.911 
12.937 
14.990 
17.068 
19.169 
21.291 
23.432 
21.193 
27.771 
29.970 

32. 181 
34.417 
36.667 

18.911 

41.216 

43.111 
41.811 

48.161 

10.112 

12.676 

11.217 
17.613 
60.066 
62.491 
64.939 
80.813 
83.163 
81.533 
87.881 
90.213 
92.SH3 

94.933 
97.281 
99.611 

101.981 

104.111 
106.681 
109.011 

111.181 
113.711 
116.001 
118.431 
120.783 
121.133 
123.463 
127.831 

110.181 

132.111 
134.861 
137.211 
139.583 

141.913 
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298. I? 
300 . 00 
400.00 
500.00 
600.00 
700.00 
900.00 
900.00 
1000.00 
f  too. 00 
1200.00 
1  300.00 
1400.00 
1  500 . 00 
1600.00 
1 700 . 00 
1300.00 
1900.00 
2000,00 
2100.00 
2200.00 
2300.00 
2400.00 
2500.00 
2600.00 
2700.00 
2800.00 
2900.00 
3000.00 
3100.00 
3200.00 
3300.00 
3400.00 
3500.00 
3600.00 
3700.00 
3800.00 
3900.00 
4000.00 
4100.00 
4200.00 
4  300.00 
4400.00 
4500.00 
4600.00 
4700.00 
4800.00 
4900,00 
5000.00 
5100.00 
5200.00 
5300.00 
5400.00 
5500.00 
5600.00 
5700.00 
5800.00 
5900.00 
6000. 00 


66.906 
66 . 906 
67.416 
68.400 
69.522 
70.663 
71.774 
72.841 
73.856 
74.819 
75.734 
76.602 
"’7.427 
78.212 
78.960 

79.676 
80.360 
81.015 
81.644 
32.248 
82. 82{> 
83.  390 
8  3.930 
84.452 
84.957 
85.445 
85.918 
86.377 
86.823 
87.255 

87.676 
88.085 
88.484 
88.872 
89.250 
39.619 
89.980 
90.331 
90.675 
91.01  1 
91.340 
91.662 
91.977 
92.285 
92.587 
92.884 
93.174 
93.459 
93.739 
94.013 
94.283 
94.548 
94.808 
95.064 
95.315 
95.563 
95.806 
96.045 
96.281 


UNCTIONS  OF  ZrOCL(t) 
56.0  KCAL/MOLE 


66.906 
66.985 
70.761 
73.818 
76. 378 
78.576 
80.501 
82.210 
83.747 
35.  1  42 
86.420 
87.598 

83.691 
89.709 
90.663 
91.560 
92.407 
93.208 
93.968 

94.692 
95. 382 
96.042 
96.674 
97.  280 
97.863 
98.424  - 
98.965 
99.486 
99. 991 
00.478 
00.950 

101.408 
01.852 
02.284 
02.703 
03.  110 
0  3.507 
03.894 
04.271 
04.638 
04.997 
05.347 
05. 6 90 
06.024 
06. 352 
06.672 
06.985 
07.292 
07.593 
07.888 
08.177 
08.461 
08.739 
09.013 
09.281 
09.545 
09.804 
10.058 
10. 309 


12.738 

12.755 
I  3.472 
I  3.902 

14.171 

14.347 
14.469 
14.555 
14.618 
14.666 
1 4 » 702 
14.731 

14. 755 
14.773 
14.789 
14. 302 
14.81  3 
14.822 
14.830 
14.837 
14.643 
1 4  •  846 
14.852 
1 *.  856 
14.860 
14. 863 
14.  866 
14.869 
14.871 
14.873 
14.875 

14.877 

14.878 

14.880 

14.881 

14.882 

14.883 

14.884 

14.885 
14.  886 

14.887 

14.888 

14.889 

14.889 

14. 890 

14.890 

14.891 

14.892 
14.892 

14.892 

14.893 
14.  893 
14.  894 

14. 894 

1 4 . 894 

14.895 

14.895 
14.  895 

14.896 


0 

.024 
1.336 
2.  709 
4.113 
5.540 
6.981 
8.432 
9.891 
11.355 
12.824 
14.295 
15.770 
17. 246 
18.724 
20.  204 
21.685 
23. 166 
24. 649 
26. 132 
27.616 
29. 101 
30.  *86 
32.071 
33.557 
35.043 
36.530 
38.016 
39.503 
40. 991 
42.478 
43.  966 
45.453 
46.  941 
48.429 
49.917 
51 . 406 
52.894 
54. 333 
55. 871 
57. 360 
58.  848 
60.337 
61 .826 

63.  315 

64.  804 
66.293 
67.  782 
69.272 
70.761 
72. 250 
73.739 
75. 229 
76.718 
78.  208 
79.697 
81.186 
82.676 
84. 166 


0 


T 

(•tt 

THKXMODYNAKIC  FUNCTIONS  OF  ZrOCL^s) 

AHf298  “  '144-°  KOU'/HOLB 

-(f-ijit/T)  • 

_ •*> _ (UUml  f> _ (Ml/il  D _ 

■i* 

ffriliwn 

298.  15 

76.175 

76.175 

16.808 

300.00 

76.175  - 

76.279 

16.835 

.031 

«00 . 00 

76.850 

31.284 

17.909 

1.774 

500.00 

78.156 

85.352 

18.520 

3.398 

600.00 

79.648 

88.764 

18.892 

5.470 

700.00 

81.164 

91.696 

19.  I  33 

7.372 

800.00 

82.644 

94.262 

19.296 

9.294 

000.00 

84.064 

96.541 

19.41 1 

11. 229 

1000.00 

'  85.416 

98.391 

19.495 

13.175 

1  100.00 

*  86.700 

100.452 

19.338 

IS. 128 

1200.00 

87.918 

102.136 

19.607 

17.086 

1  300.00 

89.074 

103.727 

19.643 

19.049 

1400.00 

90.174 

103. 184 

19.676 

21.015 

1500.00 

91.220 

106.543 

19.701 

22.984 

1 000 . 00 

92.218 

107.813 

19.721 

24.955 

1 700.00 

93  1  71 

109.01 1 

19.738 

26.928 

1800.00 

94.083 

1  10.140 

19.  752 

28. 903 

1900.00 

94.956 

111.208 

19.763 

30.878 

2000.00 

95.794 

1 12.222 

19.773 

32.855 

2100.00 

96.600 

1  1  3. 187 

19.784 

34.833 

2200.00 

97.375 

•14.107 

19.792 

36.812 

2300.00 

98.121 

1 14.987 

19.798 

38.  792 

2400.00 

98.842 

1  13.830 

19.  904 

40.772 

2500.00 

99.538 

1 16.639 

19.810 

42.753 

2600.00 

100.210 

1 17.416 

19.814 

44.734 

2700.00 

100.862 

118.164 

19.819 

46.715 

2800.00 

101.492 

1 18.984 

19.822 

48.697 

2900.00 

102.104 

119.380 

19.826 

50.680 

3000.00 

102.698 

120.232 

19.829 

52.663 

3100.00 

103.275 

120.902 

19.831 

54.646 

3200.00 

103.836 

121.532 

19.834 

56.629 

3300.00 

104.331 

122.142 

19.836 

58.612 

3400.00 

104.912 

122.735 

19.838 

60.596 

3500.00 

105.430 

123.310 

19.840 

62.580 

3600.00 

105.934 

12  3.869 

19.842 

64.564 

3700.00 

1 06.426 

1 24.412 

19.843 

66.548 

3800.00 

106.907 

124.941 

19.  843 

68.533 

3900.00 

107.376 

123.457 

19.846 

70.517 

4000.00 

107.834 

1 25. 959 

19.848 

72.502 

4100.00 

108.282 

126.<t50 

19.849 

74.487 

4200.00 

108.720 

126.928 

19.850 

76.472 

4300.00 

109.149 

127.  395 

19.851 

78.457 

4400.00 

109.569 

127.651 

19.852 

80.442 

4500. 00 

109.980 

128.297 

19.853 

82.427 

4600.00 

1 10.383 

128.734 

19.834 

84.412 

4700.00 

110.778 

129. 161 

19.834 

86.398 

4800.00 

111.166 

129.579 

19.855 

88.383 

4900.00 

1 1  1.546 

1 29.988 

19.856 

90.369 

5000.00 

H 1.918 

1 30. 389 

19.856 

92.354 

5100.00 

112.284 

1 30.783 

19.857 

94.340 

5200.00 

112.644 

1 31 . 168 

19.858 

96.326 

5300.00 

112.997 

131.546 

19.858 

98.312 

5400.00 

113.344 

1  31 .918 

19.859 

100.297 

5500 . 00 

113.685 

1  32. 282 

19.859 

102. 283 

5600.00 

114.020 

132.640 

19.860 

104.269 

5700.00 

114.350 

132.991 

19.860 

106.255 

5800.00 

114.674 

133.337 

19.860 

108.241 

5900.00 

1 14.994 

133.676 

19.861 

110.227 

6000.00 

113.308 

134.010 

19.861 

112.213 
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AH£298  "  "l0l-°  *cja/M0LE 

-(r-naw/T)  •  Cp 

-i* 

i°S _ 

■mm*. 

240.  15 

65.273 

65.277 

12.  346 

0 

:oo.oo 

65.273 

65. 349 

12.  365 

.023 

•*CO .  00 

65.760 

69.027 

1  3,  179 

1 . 30.3 

500.00 

6-6.7  30 

72,026 

1  3.695 

2.648 

'.  00 . 00 

67.820 

74.552 

14.009 

4.034 

700,00 

68.048 

76.729 

14.222 

5.  446 

300 . 00 

70.043 

78.633 

14.  369 

6.876 

?oo.c 

71.004 

80. 337 

14.474 

8.319 

1000.00 

72.006 

3 1 . 866 

14.552 

9.770 

1  100.00 

73.048 

83. 256 

14.610 

11.229 

i :oo . oo 

7  3.052 

94.529 

1 4 . 656 

12.692 

1 300.00 

74.812 

85.  704 

14.691 

14.159 

1*00.00 

75.620 

96. 793 

14. 720 

15.630 

1500.00 

76.408 

67.810 

14.  743 

17.  103 

•  600.00 

77.150 

83. 762 

14. 762 

18.578 

1700.00 

•’7.860 

69.  >.*57 

14. 778 

20.055 

1300.00 

78.530 

90.502  " 

14,791 

21.53* 

1000.00 

79.100 

01 . 302 

14.807 

23.014 

2000. 00 

70.815 

92.06  2 

14.81 3 

24.494 

2100.00 

80.415 

92.795 

14.821 

25.976 

2200.00 

80.003 

93.475 

14.829 

27.459 

2  300.00 

81.550 

94.  1  74 

14.835 

28.9*2 

2*00 . 00 

82.088 

9*4, 765 

14.840 

30.425 

2500.00 

82.607 

:■?.  371 

14.845 

71.910 

2600.00 

8  3.  109 

05. 954 

14.850 

33.394 

2700.00 

8  3.506 

9’6. 5 1  4 

14.853 

34.880 

2300.00 

84.067 

9-7.  i.15* 

14.857 

36.365 

2000.00 

84. 524 

•97.576 

14.860 

37.851 

3000 . 00 

84.967 

99 .0-0 

14.863 

39.337 

3100.00 

85.398 

98.567 

14.866 

*0.824 

3200.00 

'85.$17 

9''.  0  39 

14.969 

42.310 

3300.00 

66.2  25 

•JiC'  ^  tf  '“‘ij. 

14.870 

43.797 

3400 . 00 

6  6 .621 

■ft o , 94  (t 

14.372 

45.284 

3500 . 00 

87.008 

too. 772 

14.874 

46.772 

3600. 00 

87.  385 

100. 791 

14.676 

48.259 

3700. 00 

S' 7  7  7 

101.1  9' 3 

14.877 

49.747 

3600.00 

88! 1 12 

101.595 

14.878 

51.235 

3000.00 

8  3.46  3 

101.981 

1  *.380 

52,722 

4000.00 

88.806 

102.  358 

14. 891 

54.210 

•*100. 00 

89.141 

102.726 

14.882 

55.699 

4200. 00 

89.468 

103.084 

14.883 

57. 187 

•*300,00 

39. 739 

10  3.474 

14. 984 

58.675 

4400. 00 

00.  1 03 

103.777 

14.885 

60. 164 

4*00. 00 

90.41 1 

104.111 

14,886 

61.652 

4600. OC 

90.712 

104.  <*33 

14. 987 

63. 141 

4  700.00 

91.008 

104.758 

14. 897 

64.629 

•*  900.00 

91 . 297 

105.072 

1 4 . 888 

66. 118 

4000.00 

91.582 

1 05. 379 

14. 889 

67.607 

5000 . 00 

01.860 

105.680 

14.899 

69.096 

5100.00 

92.  1  34 

105.975 

14.390 

70.585 

5200. 00 

92.403 

1 06 . 264 

14.890 

72.074 

5  300.  00 

92.668 

106.547 

14.891 

73.563 

5400 .  On 

92.927 

106.826 

14.991 

75.052 

5500.00 

93. 182 

107.099 

14.892 

76.5*1 

5600.00 

93.433 

107.367 

14.892 

78.030 

5700.00 

93.680 

107.631 

14.893 

79.520 

5600.00 

03.923 

107.390 

14.893 

81.009 

5000 . 00 

94.162 

108.144 

14.893 

82.498 

6000. 00 

94 , 397 

109.395 

14.394 

63.998 
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THERMO0YKAMXC  FUNCTIONS  OF  ZrOF^g) 

AH®  -  -234.0  KCAL/MOLB 

.  % 

D _ (441/-1  _ 

4m 

(k Ml/Ml) 

298.15 

72.318 

72. 318 

16.142 

0 

300.00 

72.319 

72.418 

16.  173 

.030 

400.00 

72.970 

77.263 

17.450 

1.717 

500.00 

74.239 

81.244 

18, 195 

3.502 

600.00 

75.694 

84.605 

18.653 

5.346 

700.00 

77.179 

87.504 

18.951 

7.228 

800 . 00 

78.632 

90.048 

19.  153 

9.133 

900.00 

80.028 

92.31 3 

19.297 

11.056 

1000.00 

81.360 

94.352 

19.402 

12.992 

t  too. 00 

82.627 

96.205 

19.480 

14.936 

1200.00 

83.830 

97.903 

19.541 

16.887 

1300.00 

84.974 

99.469 

19.589 

18.844 

1400.00 

86.061 

100.922 

19.627 

20.  804 

1500.00 

87.098 

102.277 

1  a.  658 

22.769 

1600.00 

88.087 

103.547 

19.684 

24.736 

1700.00 

89.031 

104.740 

19.705 

26.70? 

1800.00 

89.936 

105.867 

19.723 

28.677 

1900.00 

90.803 

106.934 

19.738 

30.650 

2000.00 

91.635 

107.947 

19.751 

32.624 

2100.00 

92.435 

103.91 1 

19.762 

34.600 

2200.00 

93.205 

109.830 

19.772 

36.577 

2300.00 

93. 947 

1 10.709 

19.730 

38.554 

2400.00 

94.663 

1 1 1 .551 

19.787 

40.533 

2500.00 

95.355 

112.359 

19. 79* 

42.512 

2600.00 

96.024 

1  1  3.136 

19.8U9 

44.491 

2700.00 

96.671 

1  1  3.9e3 

19.805 

46.472 

2800.00 

97.299 

1 14.603 

19.810 

48.452 

2900.00 

97.908 

115.299 

19.814 

50.434 

3000.00 

98.499 

115.970 

19.819 

52.415 

3100.00 

99.073 

116.620 

19.921 

54.397 

3200.00 

99.631 

117.250 

19.824 

56.379 

3300.00 

100.174 

1 17.860 

19.827 

58.362 

3400.00 

100.703 

118.452 

19.830 

60.345 

3500.00 

101.218 

119.026 

19.832 

62.328 

3600.00 

101.721 

119.585 

19.934 

64.311 

3700.00 

102.21 1 

120.129 

19.836 

66.  295 

3800.00 

102.690 

120.658 

19.8  39 

68.278 

3900.00 

10’.  157 

121.173 

19.940 

70.262 

4000.00 

103.614 

121.675 

19.841 

72.246 

4100.00 

104.060 

122. 165 

19.943 

74.231 

4200.00 

104.497 

122.643 

19.944 

76.215 

4300.00 

104.924 

123.110 

19.946 

7$. 199 

4400.00 

105.343 

123.567 

19.847 

80. 184 

4500.00 

105.753 

124.013 

19. 348 

82'.  1 69 

4600.00 

106.155 

124.449 

19.949 

84. 154 

4700.00 

106.548 

124.876 

19.950 

86. 139 

4800.00 

106.935 

125.294 

19.851 

88.124 

4900.00 

107.313 

125.703 

19.852 

90.109 

5000.00 

107.685 

126. 104 

19.852 

92.094 

5100.00 

108.050 

126.497 

19.853 

94.079 

5200.00 

108.409 

126.883 

19.854 

96.065 

5300.00 

108.761 

127.261 

19.855 

98.050 

5400.00 

109.107 

127.632 

19.855 

100.035 

5500.00 

109.447 

127.996 

19.856 

102.021 

5600.00 

109.781 

128.354 

19.856 

104.007 

5700.00 

110.110 

128.706 

19.857 

105.992 

5800.00 

110.434 

129.051 

19.857 

107.978 

5900.00 

1 10.752 

129.390 

19.858 

1-09 . 964 

6000.00 

1 1  1 .066 

129.724 

19.959 

1 1 1 . 950 
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ERRATA  IN  PREVIOUS  REPORTS 


A  number  of  significant  errors  or  omissions  in  presentation 
have  been  noted  in  the  five  quarterly  reports  issued  under  this  contract. 
Corrections,  including  those  previously  given  in  the  Fourth  Quarterly 
Report,  are  as  follows: 

First  Quarterly  (15  September  1962) 

p.  11  --  The  species  "TaCljCg)"  has  been  omitted  from  the 
list  on  this  page. 

p.  16  --  Ng  should  be  defined  as  moles  of  nozzle  material  as 
condensed  phase  per  100  grams  of  equilibrium  mixture,  including  condensed 
phase . 

p.  17  and  Figures  1*8  --  the  pressure  used  in  these  calculations 
was  1000  psia. 

Second  Quarterly  (15  December  1962) 

Volume  I: 

p.  22  --  The  curve  labelled  "AI2"  should  be  labelled  '"Al",  and 
the  curve  labelled  "Al"  should  be  labelled  "A^O";  the  unlabelled  curve 
should  be  labelled  "C2". 

j  p.  66  —  4&Hf298  for  ^5(8)  should  be  listed  as  -384  kcal/mole. 

/  p.  82  --  The  phrase  at  the  top  of  the  page  should  read  "sub¬ 

tracting  the  amount  not  vaporized." 

Volumes  II  and  III: 

In  the  detailed  presentation  of  performance  results,  the  values 
of  "CF  SEA  LVL"  and  "CF  VAC"  for  the  throat  condition  ("EPSILON"  -  1.000) 
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are  Interchanged  In  every  caae  on  the  "SHIFTING  EXPANSION"  siamnary  aheet. 
For  example,  the  aecond  table  on  page  43  of  Volume  II  ahould  read: 


EPSILON  -  -  -  -  ■ 

.  CF  SEA 

LVL 

CF  VAC 

1.000 

1.236 

1.251 

2.000 

1.428 

1.459 

3.000 

1.489 

1.533 

Third  Quarterly  (15  March.  1963) 

p.  4  --  Brackets  should  surround  the  terms  following  T,  as 

C“r  '  *298  '  C"298  (I-2,8>] 

p.  6  --  The  symbol  identification  code  should  read 

o  Southern  Research 
a  Mezakl 
x  R.  K,  Kelley 

p.  33  --  The  title  for  Figure  9  should  have  the  words 
"at  1000  psia"  added. 

p.  36  --  The  title  of  the  referenced  publication  should  read 
"Thermodynamic  Properties  of  Inorganic  Substances,  V.  High 
Temperature  Heat  Contents  of  Fifteen  Refractory  Borides." 

Fourth  Quarterly  (15  June  1963) 

p.  3  --  The  first  sentence  should  include  "HBO2"  in  the  list 
of  species,  even  though  this  system  was  not  studied  in  detail. 

p.  43  --  The  final  sentence  on  this  page  is  Incorrect,  in  that 
the  additivity  rule  does  apply  if  "Ax^"  is  taken  to  be  negative,  rather 
than  arbitrarily  zero.  The  correct  value  is: 


PHILCQ 

AERONUTRONIC  DIVISION 

Fifth  Quarterly  (15  September  1963) 

p.  5  - -  The  caption  for  the  Middle  graph  hould  read  "BF2" . 
p.  14  --  The  final  phraae  on  thla  page  should  read  "the  total 
pressure  line  and  CO  line  are  almost  lndlstlngulahal  le 

p.  19  --  The  reference  cited  for  contact  .ngle  measurements 
should  be  numbered  "9". 

Section  3,  pp.  25  ff  --  Some  aaiblgulty  Is  present  in  this 
section  concerning  the  use  of  the  terms  "Denison  moc>‘l",  "simplified 
Denison  model",  and  "unmodified  Denison  model".  As  used  in  this  section, 
"simplified  Denison  model"  and  "unmodlf led  Denlaon  m  del"  are  taken  to 
be  the  same,  that  is,  the  special  case  of  the  Denlso  model  derived  in 
the  Second  Quarterly  Report  under  this  contract.  (Th.  original, more  gen¬ 
eral  model  proposed  by  Denison  is  not  limited  to  the  assumption  of 
equilibrium  at  the  wall,  and  it  does  include  allowances  for  blowing  and 
for  heat  conduction  into  the  wall.)  The  "modified  De  .Leon  model"  referred 
to  by  inference  is  one  which  was  obtained  by  semi -eng:  Irical  modification 
of  the"simpllfied  Denison  model"  to  take  into  account  blowing,  heat  loss, 
and  non-unity  Lewis  numbers;  this  nouc.  has  not  been  reported  in  detail. 
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